CHAPTER I 



INTRODUCTION. 



In this thesis we present the results of investigations of special class of Newtonian 
dynamical systems. Saying Newtonian dynamical system we assume the system 
described by second order vectorial differential equation 

= F(r,f). 

It's easy to see that such systems describe the dynamics of a mass point of unit mass 
under the action of force F according to Newton's second law. Overhead point in 
the above equation means differentiation with respect to time variable t. 

In the set of all Newtonian dynamical systems we have chosen, as a subject for 
investigation, a subclass of systems, which possess one very simple geometric prop- 
erty. This property is that one can arrange orthogonal displacement of arbitrary 
hypersurface along trajectories of the systems we are going to consider. Being less 
specific, this means that with the use of systems from this class some special geo- 
metric transformations are determined. 

The study of curves, surfaces, and their transformations has a long history (see 
§ 1 and § 2 below). Thus, our choice of the subject for investigations is not random, 
we have made it on the base of firm historical background. 

§ 1. Historical digression on geometry of curves and surfaces. 

In all times most people considered geometry as applied science. It's not surprising 
therefore, that study of many geometrical problems was initiated by every-day needs 
of human. XlX-th century was not an exception in this sense. But now, practical 
tasks, which stimulated geometers in XlX-th century, are completed and forgot 
in most. However, their ideas and methods keep the impact to current state of 
differential-geometric science. 

In spite of diversity of interests and personal tastes of researchers, starting from 
the beginning of XlX-th century the curves and surfaces became main objects of 
investigations. Development of physical science, military tasks, and economic needs 
were stimulating factors. Among others one should mention optics, ballistics, cartog- 
raphy and geodesy. For sure, some particular problems in these fields were considered 
before XlX-th century as well. Vivid example for this is the study of evolutes and 
evolvents by Huygens in 1673 (see [Huyl]). More than twenty years later Leibniz 
substantially developed ideas of Huygens. In his works we first encounter the term 
parallel curves; by means of such curves he tried to initiate systematic study of 
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caustics on wave front sets. In his works Leibniz began to use the analysis of 
infinitesimals, basics of this method was given by him in earlier paper [Lbnl] . Con- 
temporary people thought that Leibniz method of infinitesimals is difficult and 
obscure (see [Strl]). Probably for this reason his investigations on parallel curves did 
not get as popular in his epoch, as they deserved to be. From the modern point of 
view Leibniz with his works was a precursor for many researchers in XlX-th century. 
Moreover, a lot of drawings of parallel curves made by Leibniz, though paradoxical 
it sounds, can be considered as first experiments in computer geometry. In last 
few dozens of years we observe an immense development of this branch of science 
(see [Frnl]). 

The above facts say that in XVII-th century substantial break-through in study 
of curves was obtained. However, the systematic success here became possible only 
when scientific society admitted and adopted the ideas of infinitesimal analysis. 
Methods of analysis, coming to geometry, gave new point of view to broad variety 
of old problems, and gave tools for solving absolutely new problems as well. As 
historians of mathematics say (see [Dupl]), the end of XVIII-th and the beginning 
of XlX-th century is characterized by growth of works on geometry of curves and 
surfaces. As a result of such "geometric" boom mathematical science got many new 
concepts and terms. Some of them were then forgot, but others are used up to date. 

Because of great variety of papers on geometry of curves and surfaces written 
classics of science we cannot name for sure the "fathers" of differential geometry. 
Nevertheless, nowadays it is used to think that foundations of geometry of curves 
and surfaces (in its modern form) was given by Monge in his famous paper [Mngl] , 
and by Gauss in [Gssl]. Further we shall not describe history of differential geometry 
in whole, but shall pay our attention to more special its part, the theory of curves 
cmd surfaces. Expecting possible criticism in biased citation of names and results, 
we would like to say that we have no malicious intent, the only reason why this 
might happen is because we have very restricted access to historical documents. 

§ 2. Shift transformations in some classical constructions. 

In the middle of XlX-th century differential geometry has grown up to separate 
mathematical discipline. Series of monographs published at that time say in favor 
of this. Most of these books arc manuals, where wide scope of results are gathered. 
Others are devoted to special problems. Among various problems in the theory 
curves and surfaces, theory of their transformations in surrounding space took the 
proper place. Remarkable results in this direction were obtained by Bianchi. Sophus 
Lie, Gatsidakis, Blaschke, Bonnet, Backlund, Bertrand, Darboux, Tzitzeica and oth- 
ers. They considered some simple geometric transformations of curves and surfaces 
in three-dimensional Euclidean space. These investigations have enriched geome- 
try with several new notions such as Tzitzeica surfaces, affine spheres by Blaschke, 
Bianchi transformations, Backlund transformations, permutability theorem. More- 
over, these investigations gave impetus to the development of the theory of nonlin- 
ear differential equations. As appeared, such equations are naturally bound with 
geometric transformations. Most famous are Sin-Gordon equation and Tzitzeica 
equation, they are under the study up to date (see [Terl], [Boll], [Udrl].). 
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For the sake of completeness we shall consider some examples of classical geomet- 
ric constructions, which we mentioned above. 

Bonnet transformation. First systematic study of transformations of this kind 
was initiated by Leibniz in his works concerning parallel curves on the plane. How- 
ever, formalizations of Leibniz's ideas in terms of differential geometry was done by 
famous french mathematician O. Bonnet one and a half century later. His results 
are published in series of papers [Bonl], [Bon2], [Bon3], and [Bon4]. 

Bonnet considered the transformation of surfaces, which was called by him a 
normal shift. Here is the description of his construction. Denote by 



r = r(M) 



the radius-vector of a point m on some regular surface S in three-dimensional Eu- 
clidean space. Then points of parallel surface S obtained from S by the normal shift 
are determined by radius- vectors 



f (M) = r(M) + n(M) m. 

Here n(M) is unitary normal vector to initial surface, while numeric constant m is 
a parameter of construction. Characteristic features of this transformation are that 

(1) segment MM, which binds corresponding points on parallel surfaces, is orthog- 
onal to tangent planes Tm{S) and TjQi{S): 

(2) the length of this segment |MM| does not depend on the choice of point M G S, 
it is equal to m. 

If we denote by K and H Gaussian and mean curvatures of the surface S receptively, 
then same parameters of the surface S are given by the following expressions: 



l-2mH + m^K' l-2mH + m?K' 

Anther very important feature of the construction of normal shift is that it can be ap- 
plied to arbitrary smooth surface S. However, this property of normal shift, prob- 
ably, was not considered as important by Bonnet. Following traditions of Monge's 
school, he studied only those classes of surfaces, which are often encountered in 
applications. Therefore most results of Bonnet are devoted to surfaces of constant 
Gaussian curvature. The following facts were found by Bonnet, we give them here 
as they were formulated in [Bon2]. 

Theorem 2.1. Each surface of constant Gaussian curvature K = 1/m^ are as- 
sociated with two surfaces of constant mean curvature H = ±l/m, they are parallel 
to initial surface and are places at the distance ±m apart from it. 

Theorem 2.2 (inverse conjecture). Each surface with nonzero constant mean 
curvature is associated with two parallel surfaces, one of which is of constant Gauss- 
ian curvature, while other is of constant mean curvature. 
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These two theorems are among the first results in the theory of transformations 
of constant Gaussian curvature. 

Gatsidakis transformations. Next very famous transformation of surfaces was 
obtained by Gatsidakis. These transformation binds two spherical surfaces. Let S 
be an arbitrary spherical surface with Gaussian curvature K = l/a^. In isothermal 
system of local coordinates (u, v) its first and second fundamental forms are 

ds^ = E du^ + 2 F dudv + G dv'^ , 

II = L{du^ + dv'^), here L = N, M = 0. 



Formula for Gaussian curvature yields 
L _N _l 



H = ^EG-F'^. 



The Codazzi-Meinardi equation (see [Eisl], [Norl]) in this case is reduced to the 
following system of equations for coefficients of linear element of surface 



(2.1) 



dG_ _ dE_^dF^ 
du du dv ' 



dE_ _dG_ _^dF^ 
dv dv du 



Let's consider the following special solution of the system of equations (2.1): 

(2.2) F = 0, E-G = const . 
Omitting trivial case, when constant in (2.2) is zero, we can take that 

(2.3) E^a^ cosh^ O, F ^ 0, G = a? sinh^ O. 
For the coefficients of second fundamental form then we have 



(2.4) L = N = a sinh O, cosh O, 

Let's substitute (2.3) and (2.4) into Gauss equation 



M = 0. 



LN-M^ 



(2.5) 



1 

2H [du 



F dE 1 dG 



+ 



d_ 

dv 



EH dv H du 
2 dF 1 dE F dE' 



H du H dv EHdu_ 
This leads to the following differential equation for the function O: 

d'^O d^O 



(2.6) 



du'^ dv 



+ sinhO coshO = 0. 
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Theorem 2.3. Each solution of the equation (2.6) determines some spherical 
surface of Gaussian curvature K = 1/a^ . 

Proof of this theorem is quite easy. It is based on above considerations. 
Let's differentiate the relationships (2.1) with respect to u and v, then add the 
resulting expressions. This yields the following differential consequence from (2.1): 



d'^E d^E 
■ + 



Due to (2.1) Gauss equation (2.5) can be rewritten as 



(2.8) 




dE_dG_ 
du dv 



dEdG 
dv du 



+ 



d^G d^G \ 
'^'d^J 



du^ 



Let's look attentively at the Codazzi-Meinardi equations. Note that these equations 
do not change if we exchange functions E and G and simultaneously change sign of 
the function F. The same is true for the relationships (2.7) and (2.8), and for Gauss 
equation (2.5). So we have the theorem. 

Theorem 2.4. If in tree-dimensional Euclidean space linear element of spherical 
surface in isothermal system of local coordinates {u, v) has the form 

ds^ = E du^ + 2 F dudv + G dv^ , 

then in this space there exists some other spherical surface S with the same Gaussian 
curvature. Its linear element is 

ds^ = Gdu^ - 2Fdudv + Edv^. 

This correspondence between two spherical surfaces in three-dimensional Eu- 
clidean surface Bianchi named the Gatsidakis transformation (see [Bia2]). 

Bianchi transformation. Another type of surfaces, which were intensively in- 
vestigated in classical papers, are pseudospherical surfaces. In contrast to spher- 
ical surfaces, they have constant negative Gaussian curvature K = — 1/a^. Con- 
sidering pseudospherical surfaces in the same was, as we did for spherical ones, we 
get the following expressions for coefficients of linear element: 

(2.9) E = a^ cos^ co, F = 0, G = sin^ w. 

Here a; is a function of local coordinates on the surface. It satisfies the equation 

d'^uj d'^uj 



(2.10) 



du^ dv"^ 



= smw cosw. 
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An inverse conjecture is also fulfilled: each solution of the equation (2.10) determines 
some pseudospherical surface. Coefficients of first fundamental form of this surface 
are calculated according to formulas (2.9), while second fundamental form is 

(2.11) II = —a sinco cosoj {du^ — dv^). 

Theory of pseudospherical surfaces was substantially moved forward by works of 
french mathematician Ribacour and Italian geometer Bianchi. Ribacour has proved 

the following conjecture. 

Theorem 2.5. Circles of radius a on tangent plane to pseudospherical surface 
of curvature K = —1 /a^ with centers at the point of tangency are orthogonal trajec- 
tories of infinite family of surfaces of the same curvature K = —1/a^. 

Relying on this theorem, Bianchi in [Bial] proved the following one. 

Theorem 2.6. For any surface S of the curvature K = in three-dimensio- 

nal Euclidean space there exists another surface S of the same curvature K = —I / 
such that for any point M € S a point M G S is put into correspondence, and thereby 
the following conditions are fulfilled: 

(1) the distance \MM\ is constant and equal to a; 

(2) tangent planes Ts{M) and Tg{M) are perpendicular to each other; 

(3) the segment [MM] lies in the intersection of tangent 
planes Ts{M) and Tg{M). 

Theorem 2.6 describes the correspondence of two pseudospherical surfaces in 
three-dimensional Euclidean space. Darboux in [Darl] has named is a Bianchi 
transfer mat ion . 

Bianchi transformation binds two surfaces of constant negative Gaussian curva- 
ture. It's natural to ask whether can the conditions (l)-(3) be realized on the other 
surfaces. The answer to this question was given by Sophus Lie in [Licl]. He showed 
that all three conditions (l)-(3) simultaneously can be realized only on surfaces 
of constant negative curvature. 

Backlund transformations. In further development of Bianchi's construction 

the conditions (1) (3) in theorem 2.6 became more weak. Backlund in his construc- 
tion removed the condition (2). In his construction mutual arrangement of tangent 
planes Ts{M) and Tg{M) is determined by arbitrary constant angle: 

Z(Ts(M),T^(M)) =7 = const. 

Other two conditions (1) and (3) in Backlund's construction remain same as in 
Bianchi's construction. 

Investigating new transformations, Backlund showed in [Bad] , that they are also 
realized only on surfaces of constant negative curvature. But changing the angle, 
Backlund had changed some details in Bianchi's construction. As appeared. Back- 
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lund transformation binds two surfaces with constant negative curvature 



K = - 




instead of = — in Bianchi's construction. Bianchi's construction is a partic- 
ular case in Backlund's construction, it corresponds to the choice 7 = 90° 

Lie transformations. Sophus Lie introduced one more class of transformations 
for pseudospherical surfaces (see [Liel]). It is not a part of Bianchi-Backlund hier- 
archy. But it can be used to arrange a relation between above transformations of 
Bianchi and Backlund. 

Let 5* be a pseudospherical surface in three-dimensional Euclidean space. It is 
well known that at each point M G S there are two asymptotic directions, and we 
can choose asymptotic local coordinates (a, (3) determined by these two directions. 
In such coordinates linear element is given by formula 



Suppose that to = </>(«, /?) is a solution of this equation. We can check immediately 
that the function ui = (l){am, /3/m), where m is an arbitrary constant, is also the 
solution of the equation (2.12). As was shown above, each solution of the equation 
(2.10), and hence each solution of the equation (2.12) as well, determines some 
pseudospherical surface in three-dimensional Euclidean space according to formulas 
(2.9) and (2.11). 

Thus, given an arbitrary surface S of constant negative curvature, by means the 
method just described we can construct an infinite family of surfaces Sm- This 
construction is called Lie transformation. Nowadays this transformation plays 
important role in the theory of differential equations. 

Lie transformations are determined analytically, while Bianchi and Backlund 
transformations arc determined geometrically. In spite of this difference in nature of 
these transformations Sophus Lie in [Liel] managed to prove the following theorem. 

Theorem 2.7. Let be Backlund transformation binding two pseudospherical 
surfaces S and S: S = B^{S). Then we can transform S to S by means of the 
following composition of Lie and Bianchi transformations: 



Here L^ is Lie transformation with parameter 7, L^ ^ is inverse Lie transformation, 
and BgQo is Bianchi transformation. Thus B^ = L^^ o Bqqo o L^. 



ds^ = a^{da^ + 2 cos 2a; da djS + dl3'^). 



The equation (2.10) takes the following form: 



(2.12) 



dad(3 



= smijj cosuj. 



S = L-^ O Bqqo O L^ (S). 
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Darboux transformations. Backlund transformation generalizes Bianchi trans- 
formation. But its not most general construction. Further generalization was ob- 
tained by Darboux. Darboux add two parameters a and /? describing orientation 
of the segment [MM] with respect to tangent planes Ts{M) and Tg{M). Thus, in 
Darboux construction point M on initial surface S corresponds to a point M G S, 
and thereby the following conditions are fulfilled; 

(Dl) the distance |MM| is constant and equal to a; 

(D2) tangent planes Ts{M) and T§{M) intersect and form the angle 7, which is 

constant; 

(D3) the segment [MM] form the angles a and /3 with tangent planes Ts{M) and 
Tg{M) respectively, a and /3 being constants. 
Bianchi construction and Backlund construction arc special cases in Darboux con- 
struction. Here are the values of parameters corresponding to them: 

Bianchi: a = 0, /? = 0, 7 = 90°; 

Backlund: a = 0, = 0, 7 is arbitrary; 

However, some properties of general Darboux transformations are different from 
that of Bianchi and Backlund transformations. Thus, for instance, Darboux con- 
struction is realized on the surfaces, where some linear combination of Gaussian and 
mean curvatures is constant: X\K + = C = const. This construction results 
in another surface S', where some other lineal combination of Gaussian and mean 
curvatures is constant: XiK + X2H = C = const. 

§ 3. Generalization of classical construction. 

Darboux transformation terminated the series of geometric constructions, which 
are based immediately on ideas of Bianchi construction. Afterwards we observe 
the decrease in the number of papers devoted to this subject. This fact can be 
interpreted either as a slump of interest or as a period of comprehending the results 
already obtained. Second version id more probable, since later on we observe regular 
publications in Mathematical Transactions (see [Bia4]). 

Nowadays the transformations considered by Bianchi, Backlund, and Darboux 
do not loose their importance. New wave of interest to these constructions is due 
to the advent of inverse scattering method. In rather recent papers [Tenl], 
[Terl], [Chel] by Tenenblat, Terng, and Ghern we find the generalization of classical 
constructions for the case of multidimensional spaces and submanifolds in them. The 
generalization for the case of surfaces immersed into the three-dimensional space of 
constant sectional curvature was already suggested by Bianchi (see [Bia3]), this 
construction then was generalized for higher dimensions in the paper [Tenl] by 
Tenenblat. 

§ 4. Transformations of surfaces and dynamical systems. 

TO BE COMPLETED LATER 



CHAPTER II 

DYNAMICAL SYSTEMS ADMITTING 
THE NORMAL SHIFT OF CURVES IN 

§ 1. The construction of normal shift of curves 
along trajectories of Newtonian dynamical system. 

Let's consider the space equipped with the structure of Euclidean space with 
standard scalar product 



(LI) 



(X|Y) 



2 

1=1 



Suppose that 7 is an arbitrary smooth regular curve in R'^. One can assume that it 
is given in vector-parametric form 



(L2) 



vis) 



r2(s) 



where parameter s runs over some open interval of real axis, for instance, one can 
assume that s G (— 2e,+2e). It is convenient to consider closed parts of curves 
(1.2), restricting the domain of s by some segment, for instance s G [— e,+e]. In 
what follows, saying "curve" , we shall keep in mind such closed segments of smooth 
regular parametric curves: 



(1.3) 



r = r(s), s e [— e, +e]. 



On the curve (1.3) two vector-functions are defined, these are tangent vector t(s) 7^ 
and unitary normal vector n(s). 

Suppose that some Newtonian dynamical system in is given. It is defined by 
second order differential equation 



(1.4) 



f = F(r,r), 



where overhead point means differentiation with respect to parameter t being the 
time variable. The equation (1.4) describes a dynamics of mass point with unit mass 
in the force field F, according to Newton's second law. Each trajectory of dynamical 
system (1.4) is determined by fixing radius- vector r of the point at initial instant of 
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time t = and by fixing vector of velocity v = r at t = 0. Suppose that at initial 
instant of time from each point of the curve 7 a particle of unit mass is launched in 
the direction of normal vector n(s). The modulus of initial velocity for such particles 
can be given by some smooth function i' = iy{s). This is equivalent to setting up the 
following Cauchy problem for the equation (1.4): 

(1.5) r| = r(s), r| = !/(s) • n(s). 

Solution of Cauchy problem (1.5) is a vector-function 

(1.6) r = r(.s,t) 

that depend on time t and on the parameter s on initial curve 7. For each fixed t from 
some semiopen interval [0, S) the function (1.6) can be treated as a vector-function 
of one variable rt(s) = r{s,t). It determines a curve 74 in the space M^. 

The point with radius- vector rt{s) on 7^ corresponds to the point with radius- 
vector r(s) on 7 (these two points lie on the same trajectory of dynamical system 
(1.4) determined by Cauchy problem (1.5)). So we have a map ft' 7 ^ 7t, which is 
diffeomorphism for any fixed t from some sufficiently small interval [0,(5). The map 
ft'- 1 ^ 7, is a shift of curve 7 along trajectories of Newtonian dynamical 
system (1.4). 

The family of maps {/*}, where t G [0, i5), is a natural pretender for being the 
generalization of classical Bonnet construction. But to be full-value generalization 
it should preserve two properties peculiar to classical Bonnet transformations (see 
theorem 1.7 in Chapter 1). First of these properties (when transferred to the present 
situation) means that segment of trajectory binding two corresponding points A and 
Af = ft{A) on the curves 7 and 7^ should be perpendicular to these curves. 

Definition 1.1. The shift of curve 7 along trajectories of Newtonian dynam- 
ical system is called a normal shift, if all curves 7t obtained by this shift are 
perpendicular to the trajectories of shift. 

Below we consider few examples in order to understand how restrictive is the 
condition of normality formulated in definition 1.1. 



Examples illustrating the construction. 



Example 1. The shift of rectangular segment in homogeneous field of 
gravity. Suppose that we have homogeneous force field in the plane. We assume it 
being unitary and directed downward along the OY axis (as shown on Fig. 1.1): 



(1.7) 



Dynamical system (1.7) models the motion of bodies in the field of Earth's gravity 
near its surface. As an initial curve we take some part of OX axis 



(1.8) 



r{s) 



SG [-1,1], 
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Then we consider two cases. In first case wc set up Caiichy problem for the 
dynamical system (1.7) so that the modulus of initial velocity is equal to unity at 
all points of the segment (1.8): 



(1.9) 





s 


L=o 











L=o 


1 



This corresponds to the choice i/(s) = 1 in formulas (1.5). Solution of Cauchy 
problem (1.9) for the equation (1.7) has the form 



(1.10) 



r{s,t) 













s 








t + 




1 




1 








This solution describes translational motion of the chosen part of the OX axis 
downward, when it remains parallel to its initial position. Trajectories of shift 

are vertical straight lines. Therefore the normality 
condition from the definition 1.1 for the shift given 
by the formula (1.10) is fulfilled. 

In the second case we set up another Cauchy 
problem for the same dynamical system (1.7), and 
for the same part of OX axis. Let the modulus of 
initial velocity of points be the linear function of the 
parameter s in the segment (1.8) of OX axis: 



(1.11) 





s 


., 3-s 










, r = 


1 



This corresponds to the choice ^{s) = (3 — s)/4 in formulas (1.5). The solution of 
Cauchy problem (1.11) for the equation (1.7) has the form 



(1.12) 



r{s,t) = 



1 





n 4 - S 







s 










t + 




2 


1 


' - 3 


1 








Trajectories of dynamical system in this case arc vertical straight lines as before; 

however, during the shift the segment of OX axis is mapped to the segment of 

straight line such that the tangent of its slope angle 
depends linearly on t (see Fig 1.2). Therefore the 
normality condition stated in definition 1.1 is broken, 
i. e. shift (1.12) is not a normal shift. 

The above two cases in the example 1 show that 
whether the normality condition is or is not fulfilled 
for the particular shift depends drastically on a 
choice of the value of initial velocity at the points of 
curve 7 to be shifted. 
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Example 2. One-dimensional harmonic oscillator in R^. Now suppose that 
non-homogeneous force field on the plane is given, so that it is directed along OY 
axis; and suppose that its value depends linearly on the second component y = 
of radius-vector r. Then we have 



(1.13) 



r = -w y 



Dynamical system (1.13) describes oscillating motion with the frequency uj along OY 
axis and free motion along OX axis. Its trajectories in general case are sinusoids. 

In first case as an initial curve 7 we choose the part (1.8) of OX axis (just the 
same as in example 1). On this part of OX axis we set up the following Cauchy 
problem for the differential equation (1.13): 



(1.14) 





s 


lt=o 











L=o 


1 



This correspond to the choice i^(s) = 1 under the assumption that we take upper 
unit normal vector of the segment (1.8). The solution of Cauchy problem (1.14) is 



(1.15) 



r{s,t) = - 

UJ 



UJ s 

sm{u)t) 



Solution (1.15) describes oscillating motion of the segment of OX axis along another 
axis OY such that this segment remains parallel to its initial position. Trajectories 
of shift are segments of vertical straight lines (same as on the Fig. 1.1). Thus, here 
the choice of modulus of initial velocity iy{s) = 1 provides the normality condition 
from definition 1.1 to be fulfilled. 

In the second case as an initial curve 7 we choose the part of tilted line with 
slope angle 45°. It is given by 



(1.16) 



r(s) = ^ 



se [-1,1]. 



On a segment (1.16) we consider the Cauchy problem (1.5) for the equation (1.13) 
without fixing particular choice of the function i'(s): 



1 ^ 


1 


•1 ^(s) 


-1 




1 


r = — — — 
lt=o \/2 


1 



(1.17) 



The solution of Cauchy problem (1.17) then is given by 



(1.18) 



r{t,s) 



V2 



U) 



—u{s)u)t + SU) 
suj cos{ujt) + ^{s) sm{ujt) 
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Let's introduce the following function^ of the variables s and t: 
(1.19) ^it,s) = {Tt,Tt), 

This function is a scalar product of tangent vector Tt = dr{t,s)/ds to the curve 
7t at the point defined by the value of s and the vector Tt — dr{t, s)/dt tangent 
to the trajectory of dynamical system (1.13) crossing the curve 74 at that point. 
The function (1.19) characterizes the angle between the curve being shifted and 
trajectory of shift. When it vanishes, this means that the condition of normality 
from the definition 1.1 is fulfilled. In our particular case (for the shift defined by 
(1.17)) the function has the following explicit form: 

, , dvis) ( vis) dvis) \ / X . / N 

(p = vis) 7-^ t + — — 7^ - sw \ cosiLot) smiojt) + 

ds \ u> ds J 

vi^s) + s ^Mll] cos^iujt) - (vis) + s ^^^^ 
ds J \ ds 

From this formula we see that the condition of identical vanishing of ip is equivalent 
to the system of three equations for v{s): 

, , dvis) ^ vis) dvis) „ , , dvis) 

i/(s)— ^--=0, -^—^-sijj = ^, i/(s) + s — r-^ = 0. 

ds u} ds ds 

Simple analysis of this system shows that it is not compatible. Indeed, from the first 
of the above three equation we derive 

Ms) ^ Q 

ds 

Hence v{s) = const. Substituting v{s) = const into the second equation, wc get 
SLO = 0. But this cannot be valid for all s G [—1, 1], since ui ^ 0. This contradiction 
shows that on the curve (1.16) one cannot choose initial data (1.17) that provide 
normality condition stated in definition 1.1. So the shift (1.18) of the curve (1.16) 
can never be a normal shift. 

§ 2. Dynamical systems admitting the normal shift of curves. 

The examples, which wc consider in the above section, show that normality of shift 
of plane curve 7 along trajectories of a given dynamical system depends on proper 
choice of modulus of initial velocity v{s) on 7. However, from the same examples 
we see that for some dynamical systems (c. g. for harmonic oscillator) this proper 
choice is possible not for all curves. The shift along trajectories of such systems 
couldn't be full- value generalization of Bonnet construction, we loose the property 



^In thesis [Shr5] this function was called the function of deviation. Such functions play important 
role in the theory of dynamical systems admitting the normal shift. 



24 CHAPTER II. SYSTEMS ADMITTING THE NORMAL SHIFT OF CURVES. 



of universality (see theorem 1.7 in chapter I). In order to escape this situation we 
should also do the proper choice of force field of dynamical system 

(2.1) f = F(r,r). 

Thus the problem of describing the class of Newtonian dynamical systems admitting 

the normal shift arises. 

Definition 2.1. Newtonian dynamical system (2.1) is called a system admit- 
ting the normal shift in M? if for any smooth curve 7 given parametrically by 

vector- function r(,s), s G [— e, one can find a scalar fimction i'(s) 7^ such that 
the shift along the trajectories of dynamical system (2.1) determined by the solution 
of Cauchy problem with initial data 

(2.2) r| = r(s), r| = iy{s) ■ n(s) 

u=o lt=o 

is a normal shift of the curve 7 in the sense of definition 1.1. 

Definition 2.1, which was first formulated in paper [Bol3] (sec also preprint [Bol2]), 
was the starting point in constructing the theory of dynamical systems admitting 
the normal shift. Further it will be slightly modified in accordance with later results 
(see definition 2.3 below). 

Phase space of Newtoniem dynemiical system. 

Suppose that in two-dimensional space M = a Newtonian dynamical system 
(2.1) is given. The equation (2.1) can be written as a system of two equations: 



(2.3) 



F(r, 



Vector r determines the point of the space M = R^, it is called a configuration space 
of dynamical system (2.3). Pair of vectors (r, v) determines the point of the space 



(2.4) 



TM = M( 



It is called a phase space of dynamical system (2.3). Second summand F = in 
direct sum (2.4) in what follows will be called a space of velocities. The trajectories 
r = r(t) of dynamical system (2.3) and the curves r — r(s) to be shifted along trajec- 
tories of this system are assumed to be lying in the first component of the expansion 
(2.4). Tangent vectors to these curves, vector of velocity v, and vector of force F 
are referred to the second component of the expansion (2.4). The isomorphism of 
spaces M = and V = allows us to carry scalar product from M to the space 
of velocities V. Therefore in V we have scalar function 



(2.5) v=\v\ = V(v, v), 

being equal to the modulus of velocity vector, and we have vector-function 



(2.6) 



N 



V 



Ar2 
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the values of which being formed by unit vectors directed along the vector of velocity 
V. Configuration space M = is orientable, we can fix some orientation in it. This 
orientation is naturally carried to the space of velocities V = M?. Therefore we can 
consider the operator of rotation by the right angle 90° in V, and can apply this 
operator to the unit vector N in (2.6). As result we get the vector 



(2.7) 



M 





1 






~W\' 





Unitary vectors (2.6) and (2.7) form orthonormal frame with right orientation in the 
space of velocities y = R^. 

Let's consider the force field r(r, v) of dynamical system (2.3). This is vectorial 
function on a phase space, its values are given by force vector F, which is assumed 
to be an element of the space of velocities V = M^. Force vector F = F(r, v) can be 
expanded in the base formed by vectors N and M: 



(2.8) 



F = A{r, v) • N + B{r, v) • M. 



Here N = N(v) and M = M(v). It is convenient to add formal dependence of the 
vector r, considering N and M as vector-functions on phase space (2.4): N = N(r, v) 
and M = M(r, v). Therefore in the theory of Newtonian dynamical systems in 
M = we have three vector-functions N, M, and F on the phase space TM with 
the values in the space of velocities V. According to the terminology suggested in 
[Bol6] (see more details in thesis [Shr5]), such functions are called extended vector 
fields. Coefficients A and B in the expansion (2.8) are scalar functions on phase 
space TM. According to the same terminology from [Bol6], they are called extended 
scalar fields. 

Let's consider some arbitrary extended scalar field / = /(r, v). This is the func- 
tion of four variables r^, r^, v^, and v^. Partial derivatives of / form two vectors: 



(2.9) 



V/ 



dl 
dl 

Qj.2 



V/ 



First vector V/ in (2.9) is naturally called spatial gradient, second vector V/ is called 
velocity gradient of the scalar field /. Let's expand spatial and velocity gradients of 
the fields A and B from (2.8) in the base of vectors N and M: 



(2.10) 



\/A = aiN-|-a2M, 
VB = /3i N -h /32 M, 



VA = asN + a4 M, 
VB = l33N + (34 M. 



Coefficients ai, a2, as, a^, and /32, /Ss, Pa from (2.10) do characterize the force 
field F of dynamical system (2.3). If force field F(r, v) is given explicitly, then 



26 CHAPTER II. SYSTEMS ADMITTING THE NORMAL SHIFT OF CURVES. 



coefficients ai, a2, 0:3, 04, and (3i, P2, Ps, 04 also can be calculated in explicit form: 
(2.11) 



ai = (VA,N), 
a2 = (VAM), 



a3 = (VA,N), 
ai = {VA,M), 



(2.12) 



/3i = (VS,N), 
/32 = (VB,M), 



/33 = (VS,N), 
p4 = {VB,N). 



Operators of spatial and velocity gradients can be applied not only to scalar fields, 
but to vectorial fields as well. Suppose that X = X(r, v) is an extended vector field: 



X = 



Xi(r,v) 
X2(r,v) 

With partial derivatives of its components we form two matrices: 



(2.13) 



VX 





dx^ 






dX^ 


dx^ 




Qj.2 



VX 



dx^ 


dX^ 






dx^ 


dX^ 







Components of matrices (2.13) define two extended tensor fields VX and VX of the 
type (1, 1). Spatial gradients of the fields N and M from (2.6) and (2.7) are zero: 

(2.14) VN = 0, VM = 0. 

While their velocity gradients can be calculated in explicit form; 



(2.15) 



ViM^ = - 



k 



In formula (2.15) by Mi wc denote covariant components of vector field M, they are 
obtained from contravariant components of this field by lowering index: 



Mi 



Though, since metric tensor of standard Euclidean metric of in natural Cartesian 
coordinates is expressed by unitary matrix 



9ij 



1 for i= j, 
for i ^ j, 
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we can do no differences between covariant and contravariant coordinates for all 
tensor fields. Particularly, for components of vector field M we have 

Ml = M\ M2 = M^. 

Let r = r(t) be a vector-function determining some fixed trajectory of dynamical 
system (2.3). At the points of this trajectory the vector of velocity is a time derivative 
of radius- vector: v = r. Substituting r = r(t) and v = i(t) into the arguments 
of extended scalar field, we get a scalar function of t. It is easy to calculate the 
derivative of such function. When applied to the field (2.5), upon substituting 
r = v{t) and v = r{t), we get the following formulas for the derivatives: 

The relationships (2.16) follow immediately from the equation of time dynamics (2.3) 
due to the expansion (2.8). Similarly, for scalar fields A and B, upon substituting 
r = r(t) and v = r(i), we can calculate their time derivatives. Prom (2.3) and from 
the expansions (2.8) and (2.10) we get formulas 



(2.17) 



dA , , ^ „ 
— = ai |v| + as A + 0:4 B, 

^ =/3i|v|+/33A + /34B, 



where coefficients a\, as, a4, (3^, and (3i are calculated according to the formulas 
(2.11) and (2.12). 

When substituting r = r(t) and v = v{t) into the arguments of extended vector 
fields, we turn them into vector-valued functions of scalar argument t. In the case 
of vector fields N and M, one can explicitly calculate their time derivatives: 

,„ rfN EM dm BN 

^^■^^^ -dt-^' -df--^- 

The relationships (2.18) follow from (2.14), (2.15), (2.16), and from the equation 
of dynamics (2.3). In general case for the arbitrary tensorial field the following 
proposition holds. 

Proposition 2.1. For the extended^ tensor field X of the type (r, s) its time 
derivative on trajectories of dynamical system (2.3) is determined by formula 

2 2 

(2.19) = Y.^ykX-^::X+Y.P'^kX-^::X, 

fe=i fe=i 



^Recall that extended tensor field in M = is understood as multidimensional indexed array, 
components of which are the functions of double set of arguments (r, v) = (r^, r^, v^, u^). This 
concept becomes more meaningful in geometry of Riemannian and Finslerian manifolds (see thesis 
[Shr5] and Chapters III and IV below). 
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where Vk^Y'"]'' and VkX]^"'Y are components of spatial and velocity gradients of 
the field X respectively. 

Formula (2.19) expresses the rule for differentiation of composite function, since 
yk _j.k^ pk _ yk g^j^^ £qj. gp^^ial aucl velocity gradients we have 

(2.20) v.Xi-:;i: = ^i-^, VkX^lz'i:^^- 

When applied to the fields |v|, |v|~^, A, B, N, and M, from (2.19) we easily derive 
the relationships (2.16), (2.17), and (2.18). 

Vector of variation of trajectory. 

Let's consider smooth one-parametric family of trajectories^ of Newtonian dy- 
namical system (2.1) in M^. It is determined by vector- function of two variables 
r = r{t, s) that satisfies the equation (2.1) with respect to variable t. Let's calcu- 
late the following derivative 

Vector T determined by formula (2.21) is called the vector of variation on the 
trajectory rs{t) = r{t, s). Its components satisfy second order differential equations 

2 2 

(2.22) f* = ^ViF'=T*-h^ViF'='rS A; = 1,2. 

i=l i=l 

Here overhead point means differentiation with respect to t for fixed s. The equations 

(2.22) are obtained by direct differentiation in s of the equation of dynamics = 
F(r,r) written in coordinates. 

Let's consider two functions (p{t, s) and ip{t, s) obtained as scalar products of the 
above vector r with vectors N and M respectively: 

(2.23) ^ = (t, N>, ^ = (r, M). 

Functions (p(t, s) and ijj(t, s) can be determined as coordinates of vector t in or- 
thonormal right-oriented base formed by vectors N and M: t = ip-N + M. Note 
that if is an analog of the function of deviation considered in thesis [Shr5]. The 
latter is defined as scalar product (r, v). 

From (2.22) we derive differential equations for (f and ip. Let's differentiate (2.23) 
with respect to variable t and take into account the relationships (2.18): 

(2.24) (f, N) = 0-^V, ('^' M) = V; + ^ 



^An example of such family of trajectories is given by the solution of the Cauchy problem (2.2) 
set up at the points of some curve 7. 
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29 



= (f, N) + (f, N) = (f, N) + ^ ^ + ^ ^, 

(2.25) 

d(f, M) , V / ■ V / V B 

= (f, M) + (f, M) = (f, M) - - ^+ — V. 

Differentiating right hand sides of the relationships (2.24) and taking into account 
formulas (2.16) and (2.17), we get 

OT |v| |vp |v| |v| 

(2.26) 

9(t, M) .. B . ^ ^ ^ B 

=V + - 137 + PI + P3 — + P4 — 



Comparing first relationships in (2.25) and (2.26), we derive the following formula: 

Similarly, comparing second relationships in (2.25) and (2.26), we get formula for i/;: 

(2.28) v; = -^¥'+(t,M) + ^^ + (f^ -a^-p^A. 

|v| |v|2 V|v|2 |v| Iviy 

Formulas (2.27) and (2.28) contain the entry of second derivative of the vector of 
variation f. It can be expressed through t and t due to (2.22). In order to do it we 
transform the equations (2.22), which are written in coordinates, to vectorial form. 
For components of gradients VF and VF due to the expansion (2.8) we have 

^iF^ = ViA N'^ + ViBM'', 
v.f.= (v.^-^)iV'+(v.B + iiii)M^ 

In deriving (2.29) wc took into account (2.14) and (2.15). Let's multiply first equality 

(2.29) by T* and sum up in i. Thereby we use relationships (2.10): 

2 
i=l 
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Then let's multiply the equality (2.29) by f' and sum up in i. Thereby we use 
relationships (2.10) and (2.24): 



k 



+ (/33(^-^v) + (/34 + ^)(^ + ^^))M^ 



Let's add the above two relationships and collect similar terms. Now we are able to 
write vectorial form of the equation (2.22): 



V / \ V J \ V 



•V J \ V / \ V 



Written in this form, the equation (2.22) can be easily used to derive formulas for 
scalar products (r, N) and (r, M): 

(f, N) = ai ip + a2iJ + as (if - ^^'ilij + (a4 - (^+ ^^ipj, 

(f, M)=/Ji^ + /32V + /33(</j-^^) + (/?4 + ^)(V'+^<p)- 

Now substitute the above expressions into formulas (2.27) and (2.28). As a result 
we obtain two differential equations for ip and tjj. Here is first of them: 

B ■ BA A B 

<^ - TT ^ + "r~f2 V'-/3iV'-/33-r-rV'-/34T-r^ = 
|v| |vp |v| |v| 

(2.30) 

= ai if + a2il} + as [>p - -r—^ tp] +Q!4(V) + — 



Second differential equation for the functions and V' has the following form: 

B B A A B 

(2.31) ^ 
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Normal shift and differential equation 
for the modulus of initial velocity. 

Function is one of two components of the vector of variation of trajectory r 
expanded in the base of vectors N and M: 



(2.32) T = (pN + V'M. 

Differential equation (2.30) for the function Lp from (2.32) was derived in the case of 
arbitrary one-parametric family of trajectories r = r(s, t) of dynamical system (2.1). 
Now suppose that vector-function r = r(s, t) corresponds to the normal shift of some 
curve 7 along trajectories of dynamical system (2.1). In this case it is defined as the 
solution of Cauchy problem 



(2.33) r| = r(s), r| = v{s) ■ n(s), 

li=0 lt=0 

where r = r(s) is parametric equation of the curve 7, n(s) is unitary normal vector 
on 7, and v{s) is some scalar function determining modulus of initial velocity. 

For any fixed t vector-function r = (s) = r(t, s) determines in parametric form 
the curve 7* obtained as a result of shift ft - ^ Jt (see § 1 above). Therefore vector 
of variation r = T{t, s) determined by formula (2.21) is a tangent vector to curves 7t. 
The normality condition for shift (see definition 1.1) is equivalent to orthogonality 
of T and the velocity vector 

v = r=^ = |v|.N, 

which is tangent to trajectories of shift. Hence the function 

(2.34) <^ = <T,N) 

should vanish on any curve 7f . In other words, function ip should be identically zero: 
(fi{t, s) = 0. Initial data (2.33) provide vanishing of (p iov t = 0: 

(2.35) ip\ = 0. 

Let's calculate time derivative of the function (p. In order to do this we difii'erentiate 
with respect to t the expression (2.34): 

. dp d{T,N) Idr \ I dn\ 

In this expression the derivative dr/dt should be replaced by the derivative of ve- 
locity vector v{t, s) with respect to parameter s. Indeed, we have 

dr _ d / dr{t,s) \ _ d_l'dr{t,s)\ _ dw 



dt dt\ ds J ds\ dt J ds 



32 



CHAPTER II. SYSTEMS ADMITTING THE NORMAL SHIFT OF CURVES. 



The derivative dN/dt is determined from formulas (2.18). Therefore we write 

Here r = r{t, s) and v = f = v(t, s). Vectors N and M are also functions of t and 
s due to formulas (2.6) and (2.7). Let's substitute t = into (2.36). The direction 

of normal vector n on 7 and natural parameter s on 
this curve can be chosen so that for t = vector N 
coincides with normal vector n = n(s), while vector 
M coincides with unitary tangent vector r = r(s) 
on the curve 7 (sec Fig. 2.1). For vector v at initial 
instant of time t = from (2.33) we derive 

(2.37) v(0, s) = z/(s) • n(s). 

Let's differentiate the equality (2.37) with respect to 
s and take into account the Frenet equation for the derivative n'(s) (see [Norl] or 
[Shr4]): 



u'is) ■ n(s) ± k{s) ■ t(s). 



(=0 



Here k{s) is a curvature of the curve 7 at the point r = r(s). Let's substitute (2.38) 
into (2.36), taking t = thereby. As a result we get 

(2.39) (p\ = — 1 -— . 

lt=o ds z/(s) 

In the situation of normal shift from identical vanishing of ip we derive vanishing of 
its time derivative for t = 0. Due to (2.39) this leads to the following equation for 
the function i>[s) determining modulus of initial velocity: 

(2.40) ^ = B{r{s),i^{s)n{s)) ^ 

ds iy{s) 

Right hand side of the equation (2.40) depends upon the shape of curve, where 
Cauchy problem (2.33) is set up, and upon the function B = B{r,v) from the 
expansion (2.8) of the force field of dynamical system (2.3). 

Differential equation (2.40) for the function i^{s) from (2.33) is equivalent to the 
following two relationships for the function ip: 

(2.41) ifl =0, <f\ =0. 

lt=o lt=o 

The relationships (2.41) arc necessary conditions for the shift / : 7 — > 7* initiated by 
(2.33) to be normal in the sense of definition 1.1. However, these relationships (as 
well as the equation (2.40)) are not sufficient for to provide the normality of shift of 
the curve 7. 
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Weak normality condition. 

Though the conditions (2.41) are not sufficient for to provide the normality of 
shift of the curve 7, they indicate one way for building such sufficient conditions. 
Indeed, they have the form of initial data in the Cauchy problem for some ordinary 
differential equation of the second order. If for some reason the function ip{t, s) 
satisfies some linear ordinary differential equation of the second order 



(2.42) (fi = A{t) if + B{t) if 

for any fixed value of s on 7, then the relationships (2.41) provide identical vanishing 
of the function (f{t, s) on all trajectories of shift / : 7 — > 74. This fact leads to the 
following definition. 

Definition 2.2. Say that Newtonian dynamical system (2.1) in satisfies 
weak normality condition^ if on any trajectory of this system for any vector- 
function T{t) satisfying the equations (2.22) the scalar function ip = (r, N) is a 
solution of some linear homogeneous differential equation of the second order (2.42) 
with coefficients depending only on the trajectory r{t) that we have chosen. 

Let's fix some trajectory r{t) of dynamical system (2.1) and let's consider the 
equations (2.22) corresponding to this trajectory. This is the system of two linear 
homogeneous ordinary differential equations of the second order with respect to 
components of vector of variation r(t). Solutions of such system form linear vector 
space of the dimension 4. Let's denote it by T. 

Functions and tp given by scalar products (2.23) are coordinates of the vector 
of variation r in orthonormal frame formed by vectors N and M (see the expansion 
(2.32) above). They satisfy the system of two linear homogeneous differential equa- 
tions of the second order (2.30) and (2.31). The space of solution of such system is 
also 4-dimensional. It is isomorphic to 1, the expansion (2.32) written in coordinates 



(2.43) T^=(pN^ + ^M\ t"^ = ipN^ + ^M^, 

establishes such isomorphism. Hence, due to relationships (2.43), each solution of 
equations (2.30) and (2.31) determines some solution of equations (2.22). 

If weak normality condition from definition 2.2 is fulfilled, then it adds one more 
equation for the function 95. Let's express second derivative ip trough 95 and ip by 
means of (2.42) and substitute it into the equation (2.30). The result is written as 

{A - as) + - ai - a4 (p- ^4 + ^ + 

(2.44) 

'BA ^ ^ A ^ B B\ , ^ 

^ - /?! - /33 ^ - /?4 ^ - a2 + ag ^ V = 0. 

V V V / 



^Dynamical systems satisfying the condition from definition 2.2 were first considered in pa- 
per [B0I3] (see also preprint [BoI2]). However, the term "weaJc normality condition" itself was 
introduced later in paper [B0I6]. 
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Let's fix some value t = to and consider four linearly independent solutions of the 
system of equations (2.30) and (2.31). First of these four solutions we determine by 
means of the following initial data: 



(2.45) 



= 1, ifl = 0, VI = 0, VI 



= 0. 



\t=to 



\t=to 



\t=to 



\t=to 



Second is determined by another Cauchy problem with initial data 



(2.46) 



0, ^1 



0, VI 



1, i^l 



0. 



\t=to 



\t=to 



\t=to 



\t=to 



Cauchy problems for third and fourth solutions are set up similarly: 



(2.47) 
(2.48) 



1 = 0' 


(p\ 


1, 


V| = 


0, 


V 


\t=to 


\t=to 




\t=to 






1 = 0' 


lp\ = 


0, 


V| = 


0, 


V 


\t=to 


\t=to 




\t=to 







= 0, 



t=to 



= 1. 



t=to 



Due to the isomorphism (2.43) each of these four solutions of the equations (2.30) 
and (2.29) corresponds to some solution of the equations (2.22). Therefore, if weak 
normality condition is fulfilled, all these solutions satisfy the equality (2.44). If we 
substitute (2.45) into (2.44) and if wc recall that is an arbitrary parameter of 
Cauchy problem (2.45), then the equation (2.44) is reduced to 

(2.49) A -03=0. 

By the same arguments, substituting (2.46) into (2.44), we obtain 



(2.50) 



B — ai — ai 



B 



0. 



Further we substitute (2.47) and (2.48) into (2.44) and we get the equations 



(2.51) 



B 

a4 + 1— 7 = 0, 
|v| 



a2 + as ■ 



B 



0. 



Due to the equations (2.49) and (2.50) we can calculate coefficients A and B of the 
differential equation (2.42) in explicit form: 



(2.52) 



A = as, 



B = a\ + ai 



B 



The equations (2.51) were called weeik normality equations. They play more 
important role expressed by the following theorem. 
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Theorem 2.1. Newtonian dynamical system (2.3) in two-dimensional Euclidean 
space satisfies weak normality condition if and only if its force field F(r, v) sat- 
isfies weak normality equations (2.51). 

Direct proposition of theorem 2.1 was proved when we derived the equations 
(2.51). We have shown that weak normahty condition imphes the weak normahty 
equations for (f. Let's prove the inverse proposition. Suppose that the equations 
(2.51) are fulfilled. Then the equation (2.30) is reduced to 

(2.53) (fi = a3(p + (ai + a4 ^] (f. 



It's easy to see that the equation (2.53) has the same form as (2.42), coefficients 
A and B being determined by the relationships (2.52). This completes the proof of 
theorem 2.1. 

Weak normality equations (2.51) are the equations for the force field of dynamical 
system. Indeed, parameters A and B in (2.51) are expressed through F by formulas 

(2.54) A=(F, N), B = {F,M). 

This follows from (2.8). Parameters Q!2, Q^a, cti, and /^a, /34 arc determined by 
the relationships (2.11) and (2.12). Explicit form of the equations (2.51) and their 
analysis can be found in Chapter III below. 

Strong normality condition. 

Suppose that Newtonian dynamical system = F(r, f ) in satisfies weak nor- 
mality condition. Let's choose some arbitrary curve 7 determined by vector-function 
r = r(s) in natural parameter s. Then do shift this curve along trajectories of dy- 
namical system r = F(r, f ) by setting up the following Cauchy problem: 

(2.55) r| = r(s), r| = z^(s) • n(s). 

lt=o lt=o 

Regardless to the choice of u{s) in (2.55) the function ip{t, s) from (2.34) satisfies 
the differential equation of the form 

(2.56) (f = A{t) <p + B{t) if. 

This follows from weak normality condition (see definition 2.2). For the function 
u{s) we can write the differential equation (2.40): 

du{s) _ i?(r(5),K^X.)) 

This equation provides zero initial data (2.41) for the function ip{t,s): 

(2.58) (^1 =0, ip\ =0. 

u=o lt=o 
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Prom (2.56) and (2.58) wo get the identical vanishing of the function (p{t, s), which 
is equivalent to normality of shift of the curve 7. So we can provide the normality 
of shift by solving the equation (2.57). 

The equation (2.57) is an ordinary differential equation with respect to i'(s). Each 
solution of this equation is uniquely defined by fixing the value of i/(s) at some point: 

(2.59) u{s)\ = uo^O. 

\s=so 

The restriction ^-q 7^ is due to the fact that right hand side of (2.57) is singular if 
i^(s) = 0. The equality (2.59) is called the normalizing condition for the fimction 
i/(s) determining modulus of initial velocity in the construction of shift for curve 
7. According to general theory of ordinary difii'erential equations (see [Fedl]), each 

Cauchy problem of the form (2.59) for the equations (2.57) can bo solved in some 
neighborhood of the point s ~ Sq. This yields the following theorem. 

Theorem 2.2. If Newtonian dynamical system (2.1) in satisfies weak nor- 
mality condition, then for any smooth parametric curve 7 given by vector-function 

r = r(s), for any point s — sq on that curve, and for any number vq ^ one can 
mark some smaller part of this curve containing the point s = sq, and one can find 
a function ^{s) on this part normalized by the condition (2.59) such that shift de- 
termined by initial data (2.55) is the normal shift of the marked part of curve 7 

along trajectories of dynamical system (2.1). 

Let's compare the conclusion of theorem 2.2 with the definition 2.1. It's clear that 
in both cases we have the same property of Newtonian dynamical system (2.1), but 
with one remark: in theorem 2.2 we deal with some part of curve 7 instead of the 
curve itself. This remark is essential, since one can encounter the situation, when 
none of the solutions of the equation (2.57) can be extended to the whole domain 
of parameter s on 7. Such situation arises, when vector field determined by right 
hand side of the equation (2.57) is not complete on 7 (see more details in [Kobl], 
Chapter I, §1). The simplest way to exclude all such difficulties consists in slight 
modification of definition 2.1. We do it in form of the following definition 2.3. 

Definition 2.3. Newtonian dynamical system (2.1) is called dynemiical sys- 
tem admitting the normal shift in if for any smooth parametric curve 7 
determined by vector-function r = r(s) and for any point s = sq on 7 one can mark 
some part of this curve containing the point sq, and one can find nonzero function 
i/(s) on this part such that shift determined by initial data (2.55) is the normal 
shift of marked part of curve 7 along trajectories of dynamical system (2.1). 

Definition 2.3 states the condition which is a little bit weaker than the condi- 
tion from definition 2.1. It makes the construction of normal shift of curves along 
trajectories of dynamical system more local in comparison with classical construc- 
tion of Bonnet. But for our further purposes this localization is not essential. Now 

theorem 2.2 can be reformulated as follows. 

Theorem 2.3. Weak normality condition from definition 2.2 is sufficient for 
Newtonian dynamical system (2.1) inM? to admit the normal shift of curves in the 
sense of definition 2.3. 
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Theorems 2.1 and 2.3 bind the construction of normal shift along trajectories of 
dynamical system r = F(r,r) with weak normality equations (2.51) for F. However 
this bound is unidirectional: weak normality condition is sufficient, but it is not 
necessary the normality condition from definition 2.3 to be fulfilled. Just now we 
consider the example confirming this fact. It was suggested in paper [Bol7]. 

Example 3. Let f(r,v) be a force field of Newtonian dynamical system, which 
do not satisfy the equations of weak normality, for instance, this might be the force 
field of harmonic oscillator considered in example 2 (see §1 above). Let's choose 
some smooth function fj.{v) such that 



and define new force field F(r, v) = /u(|v|) • f(r, v). This new field does not satisfy 
weak normality equations as well. Therefore corresponding dynamical system does 
not satisfy weak normality condition. However, it admits normal shift of curves in 
the sense of definition 2.3. Indeed, for any curve 7 we can choose function i/(s) = 1 in 
(2.55). This means that modulus of initial velocity is equal to unity on 7, hence F = 
0. Therefore vector of velocity is constant on trajectories and |v| = const = 1. In this 
situation trajectories of shift are segments of straight lines, while the construction of 
shift in whole coincides with classical construction of Bonnet transformation. The 
normality condition for shift in this situation is obviously fulfilled. 

The idea of the example, which wc consider above, is based on the fact that some 
parts of phase space could be excluded from the process of normal shift of curves. 
The arbitrariness of curve 7 doesn't change this situation. In order to avoid such 
circumstances in [Bol7] the strong normality condition was formulated. 

Definition 2.4. Say that Newtonian dynamical system (2.1) in satisfies 
strong normality condition if for any smooth parametric curve 7 determined 
by vector- function r — r(s), for any point s = sq on 7, and for any nonzero real 
number i>q one can mark some part of this curve containing the point Sq, and one 
can find nonzero function z^(s) on this part normalized by the condition (2.59) and 
such that shift determined by initial data (2.55) is the normal shift of marked 
part of curve 7 along trajectories of dynamical system (2.1). 

The arbitrariness of curve 7 and the arbitrariness of the parameter u ^ in the 

definition 2.4 enable us to include each trajectory of dynamical system satisfying 
strong normality condition into the process of normal shift of some curve. 



H{v) = 
< i^{v) < 1 
tx{v) = 1 



for v^2, 
for 2 < u < 3 
for v^3, 



Equivalence of weak and strong normality conditions. 



Let's compare the definition 2.4 with the conclusion of theorem 2.2. It is easy to 
see that theorem 2.2 now can be formulated shortly. 
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Theorem 2.4. Weak normality condition from definition 2.2 is sufficient for 
Newtonian dynamical system (2.1) in to satisfy strong normality condition stated 
in definition 2.4- 

Theorem 2.4 can be strengthened. This was noted in paper [Shr2]. 

Theorem 2.5. Weak and strong normality conditions for Newtonian dynamical 
systems in are equivalent. 

Theorem 2.5 in paper [Shr2] was first formulated and proved for multidimensional 
case n = dimM > 2, where it states the equivalence of complete and strong nor- 
mality conditions. Peculiarity of two-dimensional case is that complete normality 
condition reduces to weak normality condition. Below we give the proof of theo- 
rem 2.5 in its specifically two-dimensional form. 

Proof. Direct proposition of theorem 2.5 follows from theorem 2.4. Let's prove 
inverse proposition. Suppose that Newtonian dynamical system satisfies strong nor- 
mality condition. Let's choose and fix some arbitrary trajectory r = r{t) of this 
system and let's mark some point po on it. Without loss of generality we can assume 
that marked point corresponds to t = 0, its radius- vector is fq = r(0). Moreover, 
we shall assume that po is regular point of trajectory r = r(t). This means that 
velocity vector vq = r(0) at the point po is not zero. Denote i^o = |vo|- And finally, 
we denote by no the unitary vector directed along the vector of velocity vq at the 
point Po on the chosen trajectory: 

(2.60) no = N(0) = ^ 

|vo| 

Suppose that 7 is some arbitrary smooth curve in R^, satisfying the following 
two conditions: 1) it passes through the marked point po on the trajectory r = r{t) 

that we fixed above; 2) it is perpendicular to the 
trajectory r = r{t) at this point. Curve 7 can be 
determined by vector-function r = r(s) in natural 
parametrization; without loss of generality we can 
assume that the point po corresponds to s = 0. Vec- 
tor no in (2.60) is the unitary normal vector to 7 at 
the point po- It determines normal vector n = n(s) 
on the whole curve such that no = n(0). 

Remembering that dynamical system r = F(r, r) 
satisfies strong normality condition, we choose the normalization for i'(s): 

(2.61) u{s)\ = I/O. 

ls=0 

According to definition 2.4, we can find the function iy(s) normalized by the condition 
(2.61) at the point po, and such that it defines normal shift of some part of the curve 
7 along trajectories of dynamical system being considered. The trajectory r = r(t), 
which we have fixed, is one of the trajectories of such shift: 

r(t)=r(t,s)| . 

ls=0 
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Apart from this trajectory, vector-function r = r{t, ,s) obtained as the solution of 
Cauchy problem (2.55) for the equation f = F(r,r) describes the whole family of 
trajectories coming out from various points of the curve 7. This family of trajectories 
determines the vector of variation 



ds 

whose components satisfy the equations (2.22). Formulas (2.23) define two functions 

tf and tp that satisfy the equations (2.30) and (2.31). Due to normality of shift first 
of them is zero identically: ip = (r, N) = 0. Therefore its derivatives are zero too: 

(2.62) ip\ =0, <p| =0, (p\ =0. 

lt=o lt=o lt=o 

First of the equalities (2.62) is fulfilled only due to special form of initial data (2.55). 
Second equality is fulfilled due to (2.39) and the equation (2.40), which determines 
the function v{s) on 7. Let's write the third equality (2.62). In order to do it we 
use the equation (2.30) and express second derivative (p: 

... { B\ ., (BA ^ 

<fi\ = a4 + — •■^ _ 



lt=o \ |v|/ U=o 
(2.63) 

A B B\ , 

|v| |v| \\\J lt=o 

Now let's calculate the values of functions tp and tj) ior t = 0. Natural parameter s 
on the curve 7 can be chosen such that for t = the vector N coincides with normal 
vector n. This follows from formula (2.60) and from the way how the normal vector 
n(s) on 7 was chosen. Therefore we get 

(2.64) V| =(t, M)| =(t, r) = l. 

Similarly we calculate time derivative of the function ip. In order to do it we differ- 
entiate in t the expression -0 = (r, M), being the definition of tp: 

It this expression partial derivative dr/dt can be replaced by the derivative of v{t, s) 
with respect to s. Indeed, we have 

dr _ d / dr{t,s) \ d / 9r(t,s) \ _ 5v 
'dt ~ dt\ ds ) ~ ds\ dt ) ~ 'ds' 
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The derivative dM/dt is determined by formulas (2.18). Therefore we write 
(2.65) ^ = (|,M)-{r,N>^. 

Second summand in formula (2.65) vanishes for t = 0. In order to transform first 
summand we use formula (2.38). This yields 



(2.66) V =(^,M 



it=o 



\ds' 



= ±k{s). 



Now we substitute (2.64) and (2.66) into the equation (2.63). Then remember that 
we have an infinite variety of curves passing through the marked point po and being 
perpendicular to the fixed trajectory r = r{t). Among them one can find a curve 
with any preassigned value of curvature = A;(0) > at the point po- Therefore the 
equation (2.63) breaks into two separate equations: 

B ^ 
a4 + — = 0, 
|v| 

(2.67) 

BA ^ ^ A ^ B B 

— - /5l - /53 1-7 - P4 -p-T - Q!2 + Q!3 1— r = 0, 

vr V V V 



These equations (2.67) are exactly the same as weak normality equations (2.51). 
Since po is an arbitrary point of M^, the equations (2.67) hold at regular points of 
any trajectory. This means that they are fulfilled at all points of phase space, where 
|v| ^ 0. According to theorem 2.1 weak normality equations are equivalent to weak 
normality condition from definition 2.2. Theorem 2.5 is proved. □ 



CHAPTER III 



ANALYSIS OF NORMALITY EQUATIONS IN 



§ 1. Transforming the normality equations to Cartesian coordinates. 

The concept of Newtonian dynamical system admitting the normal shift of curves 
in , which arose in generahzing classical Bonnet construction, leads to equations 

B 

(1.1) 

BA A B B 

|v|2 |v| |v| |v| 

They were called the weak normality equations. In two-dimensional case they are 
simply called the normality equations, since they exhaust complete^ system of nor- 
mality equations in this case. Normality equations (1.1) are partial differential equa- 
tions with respect to components of force field F of Newtonian dynamical system 

(1.2) r = v, v = F(r,v). 

However, looking at the equations (1.1), this is not so obvious. Force field F of dy- 
namical system (1.2) are represented in the equations (1.1) only through coefhcients 
A and B of the expansion 

(1.3) F = A-N + B-M, 

and through parameters a2, 03, 04, /3i, /?3, /?4, which are defined by the expansions 
of gradients VA, VB, VA, and VB: 

= ai N + a2 M, = as N -t- M, 

(1.4) 

VB ^ iSiN + VB = PsN + P4:M. 

Our nearest goal is to transform the equations (1.1) so that components of force 
field F would be represented in explicit form. From the expansions (1.3) and (1.4) 



^In multidimensional case n ^ 3 complete system of normality equations includes weak normal- 
ity equations and additional normality equations (see more details in [Shr5]). 
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we derive the following expressions for B and 04 : 

2 

(1.5) B = {F,M) = Y,M'Fi, 

i=l 

2 

(1.6) a4 = (VA, M) = ^M' ViA. 



For coefficient A in expansion (1.3) we have the expression similar to (1.5): 

2 

(1.7) A={F,N)=J2^'Fj- 

i=i 

Due to the relationships (1.5), (1.6), and (1.7) the first normality equation equation 
in (1.1) is rewritten as 

(1.8) ^ ( |v|-i + ^ Vi {N^ Fj)] = 0. 

i=l \ j=l / 

Now let's do similar transformations in each term of the second normality equation 
in (1.1). From (1.5) and (1.7) we derive 

2 2 

I I i=ij=i 

For second and fifth terms of the second equation (1.1) we have 

-/3i -a2 = -{VB, N> - (V^, M). 

This follows from (1.4) (see relationships (2.11) and (2.12) in Chapter II). Taking 
into account (1.5) and (1.7) we derive 

2 2 

(1.10) -Pi -a2 = - X! + ^^^Pj) 

i=i j=i 

Remember that ViN^ = and V jM^ = (see relationships (2.14) in Chapter II). 
We took into account these relationships in deriving (1.10). 

Now let's transform fourth term in the second normality equation (1.1). For the 
coefficient ^4 in this term we have 

2 2 

(1.11) /?4 = (VB, M) = ^ ^ V, (f, MO M^. 
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In expanding brackets under the signs of summation in right hand side of (1.11) we 
take into account the relationships (2.15) from Chapter II. This yields 

' ' j=\ r=l ' ' i=l j=\ r=l ' ' 

Further calculations are similar to those we have already done. These calculations 
do not require special comments: 



(1-13) -/33a = -EEE 



V V 



(1-14) "3kh =EEE- ■ 



|V| V 



Let's add the equalities (1.9), (1.10), (1.12), (1.13), and (1.14). Collecting similar 
terms, we take into account the relationship 



1 for j = r, 
for j ^ r. 



which follows from the fact that vectors N and M form orthonormal frame in stan- 
dard euclidean metric of the space M^. As a result we have 

2 2 

i=i j=i 

(1.15) 

i=l j=l \r=l I ' ' ' / 

Let's denote v = |v|, let's change sign in (1.15) and join the equations (1.8) and 

(1.15) into one system. This yields the system of two partial differential equations 
of the first order. They are written as follows: 

(1.16) ^ L-ij; + ^Vi(7v^>,) J M* = 0, 

i=l \ j=l J 

2 2 

E E i^i^i + ^^^i - 2 Pi Pj) + 

i=l j=l 



(1.17) 



EEf^-E^^-^l- = o. 
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These are the normahty equations for the force field F(r, v) of Newtonian dynamical 
system = F(r, f) in R^, written in Cartesian coordinates. 

In (1.16) and (1-17) we have the entries of components of unitary vector N. 
It is uniquely determined by velocity vector: N = |v|~^ • v. Unitary vector M 
is orthogonal to N, it is also determined by velocity vector (see formula (2.7) in 
Chapter II). However, in order to determine it we should fix some orientation in 
the space ffi.^. Change of orientation leads to the transformation M to — M. This 
transformation changes sign of left hand side of the equations (1.16) and (1.17), but, 
in essential, this do not change the equations themselves, since their right hand side 
is equal to zero. We can make the normality equations (1.16) and (1-17) completely 
invariant with respect to transformation M to — M. In order to do it we define the 
operator field P with the following components: 

(1.18) P; = S'j - N' Nj = M* Mj. 

Now, multiplying the equations (1.16) and (1.17) by M^, we write 



(1.19) Fi + Y.^i Fi) ]Pl = Q, 

2 2 

E E ^"^^^i + ^i^^ - 2 Fi Pi + 

i=i i=i 



(1.20) 



i=l j=l \ r=l / 



Operator P with components (1.18) is an operator of orthogonal projection to the 
line perpendicular to velocity vector v. It is invariant with respect to the change of 
vector M by -M. 

It's worth to note that the normality equations written as (1.19) and (1.20) co- 
incide by their form with the weak normality equations for Newtonian dynamical 
system on an arbitrary two-dimensional Riemannian manifold M (sec [Bol6] and 
thesis [Shr5] ) . The difference of and the case of arbitrary Riemannian manifold 
M is concentrated in the definition of differentiations V and V. This fact will be 
used below in Chapter IV). 



§ 2. Scalar ansatz. 

Let's consider first normality equation in the system (1.1). Using this equation, 
we can express coefficient B in the expansion (1.3) through a4: 



(2.1) 



B = — Ivl Ui- 



§3. POLAR COORDINATES IN THE SPACE OF VELOCITIES. 



45 



Parameter 04 in (2.1), in turn, is determined by coefBcient A in the expansion (1.3) 
according to the formula (1.6). Therefore we have 

2 

(2.2) B = -\v\J2M'S/iA. 

Substituting (2.2) into the expansion (1.3) we get formula that expresses components 
of force vector F through coefficient A: 



(2.3) Fk = ANk-\v\YM'MkViA. 

i=l 

If we take into account (1.18), then formula (2.3) can be transformed to the form 

2 

(2.4) Ffe = A7Vfc-|v| ^P|ViA 

i=l 

Formula (2.4) was called scalar ansatz. In [Bol9] (see also thesis [Shr5]) such 
formula was used for simplifying normality equations in multidimensional case n ^ 3. 
Here in two-dimensional case this formula allows us to reduce system of normality 
equations (1.19) and (1.20) to one differential equation with respect to the function 
A(r,v): 

2 / 22 

^ Vi A + I V| ^ ^ P?"^ V, A \/r\/^A - 
i=l \ q=l r=l 



(2.5) 



2 2 \ 

^ iV A Vi^ - |v| J2 Vi A 

r=l r=l / 



0. 



Here P''' = AP" are components of projection operator P. The equation (2.5) 
will be called reduced normality equation. It is as a result of substituting (2.4) 
into the equation (1.20). Substituting (2.4) into the equation (1.19) we get the 
identity. There is no need to reproduce these calculations here, since they do not 
differ from those in multidimensional case (see [Bol9] and thesis [Shr5]). 

§ 3. Polar coordinates in the space of velocities. 

In reduced normality equation we have the entries of vector N directed along the 
vector of velocity v, and wc have the entries of vector M perpendicular to v. This 
can indicate the presence of some rotational symmetry in the space of velocities. 
In order to use this symmetry for the purposes of further simplifying the normality 
equation (2.5) let's introduce polar coordinates in the space of velocities. In place 
of Cartesian coordinates and v"^ of velocity vector v we consider radial variable 
V = |v| and angular variable 6. As 6 we can choose the angle between velocity 
vector and some fixed direction in M^, for instance, with the direction positive OX 
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semiaxis. However, we consider more flexible construction, when is referenced to 
some unit vector m = m(r). Let m = m(r) be unit vector perpendicular to m: 



(3.1) 



TO 



TO 



TO 



TO 
1 



Components of m and rh depend on spatial variables 
and r^, but they do not depend on components of 
V. Hence their velocity gradients are zero: 



(3.2) 



V,TO^ = 0, 



The derivative of any unit vector is perpendicular to 
this vector. This yields the following relationships: 

(3.3) V im-' = u)i m-' , V i'ih-' = —Wi m-' . 

Functions wi(r) and a;2(r) in (3.3) are bound with each other by the relationship 



(3.4) 



Let's recalculate the reduced normality equation (2.5) to spherical coordinates in 
velocity space. This means that we should do the following change of variables in it: 



(3.6) 



= J/, 

= V (m} cos 9 + fh^ sinfl), 
v"^ = V {m? cos9 + TO^ sm9). 



Prom (3.6) we obtain the following expressions for vectors N and M: 
(3.7) 

Now let's calculate the derivatives of v in a; and y due to (3.6): 



N = m • cos ^ + ih • sin ^, 
M = — m • sin ^ + m • cos 6. 



d_ 
dx 

d_ 
dy 



= VLOi COS ( 



VUJ2 COS 9 



TO^ 



TO^ 



— VLOi sm ( 



V0J2 sin 9 











= VUJi 




vn? 




TO^ 




Ml 




= VU)2 






M2 
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These formulas are immediate consequences of (3.3) and (3.7). The derivatives of v 
in V and 9 are calculated in a similar way: 



d 




















= COS 9 




+ sin^ 








dv 










iV2 















Ml 




= V COS 9 




— V sm.9 




= V 


M2 



d_ 
89 



Now we are able to write the Jacoby matrix for the change of variables (3.6): 







(3.8) 



J = 







In order to recalculate differential operators and Vj to variables x, y, v, 9 we 
should find the matrix inverse to matrix J in (3.8): 



(3.9) 



TVi N2 

-UJl -UJ2 v~'^ Ml M2 
Using formulas (3.9) for spatial gradients, we obtain 



1 
1 



(3.10) 



„ d d 



V - — - — 



In a similar way for velocity gradients Vi and V2 we obtain 



(3.11) 



Vi = iVi — + ^ — 
ov V oO 



~ ^8M2d 

V2 = A^2 1 — • 

dv V 89 



48 III. ANALYSIS OF NORMALITY EQUATIONS IN R'^. 

Now we can do the change of variables (3.6) in the equation (2.5). First of all, let's 
rewrite this equation itself 



(3.12) 



2 / 2 2 

^ V.A + |v| 53 53 M« \/gA VrViA - 

i=l \ q=l r=l 

2 2 \ 

- 53Af''AVrV,A- \v\^N''Wr^iA \ M' = 0, 

r=l r=l / 



applying the fact that = Mfc and Pi'^ = M" M'^. Prom the relationships 

(3.11) , as an immediate consequence, we have 

Using the relationships (3.13) and taking into account formula (3.7) for the vector 
M, we calculate third summand in left hand side of the equation (3.12): 

(3.14) -tt^'^ V.V.AM^ = -A ^( ^) = + ^. 

r— 1 r— 1 ^ ^ 

Then consider second summand in left hand side of the equation (3.12): 

|v| t V,A ±±M^Mr V. V. A = ^ | (t " V ^ ^ V.^. 

q=l i=l r=l \i=l / i=l 

Upon calculating dM^/dO on the base of (3.7), for this summand we get 

2 2 2 

(3.15) \y,\Y^M'^V,A^^M'M^\7rViA=^^ + ^^. 

q=l 2=1 r=l 

Now we are to transform first and fourth summands in reduced normality equation 

(3.12) . First of all, note that 



53^'" VrM' = 0. 



r=l 



This follows from the relationships (2.14) in Chapter II. Now for the fourth summand 
in reduced normality equation (3.12) we have 



-|v| 53 53 Vr-Vi^l = -t; 53 ^^'^ V J - 

i=l r=l r=l ^ 

(3.16) 

= - ATI Ae^ - Aey + {N^ ui + L02) Aee- 
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In deriving formula (3.16) the relationships (3.10) were used. We use them for to 
transform first summand in (3.12) as well: 

2 

(3.17) '^'^ = M^A^ + Ay - (M^ wi + ws) Ae. 

i=l 

Let's add the relationships (3.14), (3.15), (3.16), (3.17) and take into account the 
relationships (3.7) for vectors N and M written in coordinates: 

cos^ [rh^{Ax — wi Ao) + m^[Ay — L02 A$) — 

- m^(Aex - oji Aee) - m'^{A0y - 012 A00)) - 

(3.18) - sin 9 {m^ (A^-LJiAe) +171^ {Ay -iO2A0) + 

+ m''^{A0x-uiiA00)+rh'^{A0y-u!2A00)) - 

AAe^ , AA0 , A0A00 A0A^ ^ 



As a result we have transformed the normality equations (2.5) to variables a;, y, w, 

d. Last two variables v and 9 are polar coordinates in the space of velocities. First 
two variables x and y are Cartesian coordinates in configuration space M = R^. 



§ 4. Some simplest solutions of normality equation. 

According to the results of Chapter II, constructing the Newtonian dynamical 
systems in admitting the normal shift reduces to solving normality equations 

(1.1) . Latter ones are reduced to the equation (2.5) for the function A. Upon change 
of variables (3.6) the equation (2.5) turns to (3.18). Formula (2.2) in variables a;, y, 

6 has the following very simple form: 

(") --^ 

(see relationships (3.13) above). For the force field F of Newtonian dynamical system 
from (4.1) and from the expansion (1.3) we derive 

dA 

(4.2) F = A • N - — • M 

06 



Simplest solution of the equation (3.18) is given by an identically zero fimction 
A{x,y.v,9) = 0. Such solution corresponds to Newtonian dynamical system with 
zero force field F = 0; normal shift along trajectories of such system coincides with 
classical Bonnet transformation. 
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Next (by complexity) solution of the equation (3.18) is determined by arbitrary 
smooth^ function of one variable A = A{v). If we substitute A = A{v) into (4.2), 
we obtain F = A{v) ■ N. Force vector is directed along the vector of velocity, its 
modulus depends only on modulus of velocity vector |v|: 



(4.3) 



F = 



^(|v|) 



Force field (4.3) should also be considered as trivial one. Trajectories of correspond- 
ing dynamical system are straight lines, normal shift along them coincides with 
classical Bonnet transformation. 

Less trivial examples of dynamical systems admitting the normal shift were con- 
structed in paper [Bol3] (see also preprintBol2). 

Example 1. Spatially homogeneous, but not isotropic force field. Let's 
choose constant vector field m in (3.1): 



(4.4) 



Then u)i and u)2 in formulas (3.3) are equal to zero. Such choice of m and m 
substantially simplifies the normality equation (3.18): 



(4.5) 



{Ay - Ae^) cose - {A^ + Agy) sme + 

AgA^ AAg 



AAg Ag Agg 



Among all solutions of the equation (4.5) now we choose those, which do not depend 
on X and y, i. e. ^ = A{v,&). This provides spatial homogeneity. Then 



(4.6) 



AAg Ag Agg Ag Ay 
o ~r o ~r 



A Agy 



Let's take A = A[v) cosO, where A{v) is an arbitrary smooth function of one vari- 
able. It's easy to check that the function A = A(v) cos 6 satisfies the equation (4.6). 
Substituting this function into the formula (4.2), we obtain 

(4.7) F = A{v) ■ (cos 61 • N + sin 6* • M). 
In Cartesian coordinates force field (4.7) is written as 

(4.8) F = A(|v|)- (2(N, m)-N-m), 



^The brecik of smoothness at the point = is admissible, since the equation (3.18) is written 
only for those points, where v ^ 0. 
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where m is constant vector, which is unitary in our example and directed along OX 
axis (see formulas (4.4)). 

As we see on Fig. 4.1, force vector in force field 
(4.7) form the angle with velocity vector v and it 
forms the angle 2 9 with OX axis. Newtonian dynam- 
ical system with force field (4.7) admits the normal 
shift of curves in the space R^. Moreover, it is in- 
tegrable in quadratures. In order to show this let's 
find its trajectories. Force F in (4.7) does not de- 
pend on coordinates and r^. Therefore the equa- 
tion = F(r, r) can be written as v = F(v). Let's 
integrate this equation in variables v and 6. Using 
the relationships (4.4), we write 





COS 9 




COS 6* 




— sin 9 


V = v ■ 


sin 9 






+ v9- 








sin^ 




cos 9 



As for the formula (4.7) for force vector, we rewrite this formula as follows: 



(4.10) 



F = A{v) cose- 



cos ( 
sinl 



+ A{v) sin9 ■ 



— smf 
cos^ 



Comparing (4.9) and (4.10), we see that vectorial equation v = F(v) is written as a 

system of two scalar equations 



(4.11) 



ii = A{v) cos 6, 
v9 = A{v) sin0. 



Solution of the system of differential equations (4.11) is the pair of functions v{t) 
and 9{t). If A{v) ^ and sin 6* 0, then function 9 = 6{t) is invertible. The 
inverse function is t = t{6). Substituting t = t{9) into v{t), we get the function 
V = v{t{9)) = v{9). From the system of two equations (4.11) we get the following 
differential equation for such function: 



(4.12) 



dv 
d0 



V A{v) cost 



= V cotan 9. 



Q v~'^ A{v) sin9 
Differential equation (4.12) is easily integrable: 
(4.13) v{9) = C I sin9\, where C = const > 0. 



Upon substituting (4.13) into the second equation (4.11) we obtain differential equa- 
tion of the first order for the function 9 = 9{t): 



(4.14) 
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The equation (4.14) is integrable in quadratures: 

^ C|sin6'|d6l 



(4.15) 



>l(C|sin6'|) sin 6* 



= t + C. 



Hence the whole system of equations (4.11) for v{t) and 6{t) is integrable in quadra- 
tures. Its solution contains two constants of integration C and C. 

Further, in order to find trajectories one should integrate the equation r = v. 
Right hand side of this vectorial equation in variables x, y, v, 6 is determined by 
formula (4.9) for the velocity vector: 



d 


X 




COS 9 






= V • 




di 


y 


sin^ 



Knowing v{t) and 6{t), we find x{t) and y{t) simply by integration: 

t t 
(4.16) x{t) =xo + jv{t) cos{e{t)) dt, y{t) = ?/o + jv{t) sm{e{t)) dt. 



Function A{v) in formula (4.7) is arbitrary one. If we choose it being constant 
A{v) = ^0 = const, then integral in (4.15) can be found in explicit form: 



(4.17) 



CO 

—— sign(sin 6) =t + C. 



Let's set up the following Cauchy problem for the equations (4.11): 

(4.18) vit)\ = vo, 9{t)\ 

\t=o 



Suppose that Aq > and > 0. Then we have 

e{t) ^Oo + ut, 



\t=0 



(4.19) 



v{t) = — sin(^o +i^t). 



Solution of Cauchy problem (4.18) is determined by formulas (4.19) only within 
some interval of values of parameter t: 



(4.20) 



TT - 6»o 

-- ^ -, 



angular velocity u) in (4.19) and (4.20) being determined by formula 
(4.21) a;= — sin6lo. 
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If we substitute (4.19) into formulas (4.16), we can calculate integrals (4.16) in 
explicit form. For x{t) and y{t) this yields 



(4.22) 



x{t)=xo-j^ (cos(2 eo + 2(jt)- cos(2 Oo)) , 
Aot Ao 



y{t) =yo + ^ - ir^ (sin(2^o + 2a;i) 



sm 



Prom formulas (4.22) we sec that if A{v) = Aq = const, trajectories of Newtonian 
dynamical system with force field (4.7) are cycloids. Thc^sc are curves drown by a 
point of a circle rolling along the straight line parallel to OY axis. 

Example 2. Spatially non- homogeneous force field with marked point. 

Lot's mark the point coinciding with the origin of Cartesian coordinates and let's 
direct the vector field m in (3.1) along the radius-vector: 



(4.23) 



m : 



1 



-y 



Here |r| = ^/x^ + y^. Using explicit formula (4.23) for vectors m and lii, we can 
calculate parameters oji and 0J2 in formulas (3.3): 



(4.24) 



m., 
|r 



^, i = h2. 



Further we use normality equation (3.18). When we take into account (4.23) and 
(4.24), this equation is written as 

{xAy-yA^-x Aqx - y Aoy - Aq) cos 6 



(4.25) 



{xAx + yAy-y Ag^ + x Aoy - Aee) sin 6 



AAqv ^ AAg ^ ApAee ^ ApA^ _^ 

With the aim of further simplification of the equation (4.25) we use polar coordinates 
in configuration space M = M^. Let's do the following change of variables: 



(4.26) 



( x = p cos 7, 
\ y = p sin7. 



Polar coordinates v and 6 in the space of velocities are not changed thereby. From 
(4.26) we derive formulas for derivatives: 



d 



(A ').7\ 



— = COS7 ■ 
OX op 



d sin 7 d 



p 97' 
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Taking into account the relationships (4.26) and (4.27), we can write (4.25) as 



{A^ - Ae- p Agp ) C0S6> Aq Agg Ag Ay 



(4.28) 



{p A I, + - . Uo) sill \o 



V 

AAe 
t;2 



If we restrict ourselves with functions A = A{p,v,6), which do not depend on 7, 
then we easily find the following solution of the equation (4.28): 



(4.29) 



A{v) cos( 



Here A(v) is an arbitrary smooth function of one variable. Applying (4.2), we see 
that the solution (4.29) corresponds to Newtonian dynamical system with force field 



(4.30) 



A{v) ■ (cos 6* ■ N ^ sill ■ M) 



In Cartesian coordinates in force field (4.30) takes the following form: 

2(N, r)-N-r 



(4.31) 



F = A(|v| 



Fig. 4.2 illustrate the geometry of force field (4.31). It is similar to geometry of the 
field (4.7). Vector of force F forms the angle 6 with vector N, and it forms the angle 
2 9 with radius- vector r, which, in turn, forms the angle 7 with OX axis. Newtonian 
dynamical system with force field (4.30) admits normal shift of curves in M^. Similar 

to dynamical system considered in example 1 above, 
this system is integrable in quadratures. In order to 
integrate it let's write the equations of dynamics (1.2) 
in variables p, 7, v, 9. For components of radius 
vector r we have the equalities 



(4.32) 



T = p- 



P- 



COS 7 
sin 7 
COS 7 
sin 7 



+ pi- 



— sin7 
cos 7 



Let's write analogous equality for components of ve- 
locity vector, using Fig 4.2 for this purpose: 





cos(7 + 9) 




cos 7 




— sin 7 


V = i; • 


sin (7 + 6) 


= V cos 9 ■ 




+ v sm9 ■ 








sin 7 




cos 7 
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Comparing this formula with (4.32), from r = v we derive two equations: 

( p = v cos 6, 
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(4.33) 



pj = v sinO. 



For the derivative of velocity vector we have the expansion, which follows from the 
fact that vector v forms the angle 'y + 6 with OX axis: 



cos(7 + 6) 
sin(7 + 0) 



- sin(7 + 6) 
cos(7 + 9) 



Force vector F in (4.30) forms the angle 29 + j = 9 + {-j + 9) with OX axis (see 
Fig. 4.2), its modulus is equal to A{v) / p. Using elementary identities from trigonom- 
etry, for the vector F we obtain the expansion 



A{v) 



cosf 



cos(7 + 9) 
sin(7 + 6) 



A{v) . 
H — sm( 



— sin(7 + 9) 
cos(7 + 9) 



In essential, it coincides with (4.30). Now from v = F we derive two more differential 
equations in addition to the equation (4.33): 



(4.34) 



A(v) 
V = cos 9, 



v{j + 9) = ^ sin^. 



Let's express 7 through p, v, and 9 by means of (4.33) and substitute the obtained 
expression into (4.34). As a result the equations (4.34) turn to 



(4.35) 



V = cos ( 



v9 



A{v) 



smf 



Similar to (4.11), the equations (4.35) can be treated as differential equations deter- 
mining function v = v{9) in parametric form: 



(4.36) 



dv V 
d9 



_ vA{v) 
Q A{v) - 



cotan^. 



Differential equation (4.36) can be integrated in quadratures: 

V 

A(v) - 



(4.37) 



vA{v) 



dv = ln(| sin^l) + const . 
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Having determined the function v{9) from (4.37), from (4.33) we derive the differ- 
ential equation for the function p{9): 

(4.38) 'la = - = TT^. -ot.ne. 

d9 A{v) - ti2 

The equation (4.38), which we have just obtained, is also integrable in quadratures: 



Let's determine the function p{6) from (4.39) and then substitute it into second 
equation (4.35) together with the function v{6) obtained from (4.37): 

(4.40) ^ = ^M:i^sin0. 

pv 

The solution of (4.40) is the function 6 = 6{t), which is calculated in quadratures 
according th the following formula: 



(4.41) I 



pie) vie) de 

> ^ > =t + const. 



A{v{e)) - v{e)'^ sine 



Having determined 9 = 9(t) from (4.41) and having substituted it into v{0) and p(9), 
we get the functions v = v{t) and p = p{t). We are to determine the last of four 
functions v{t), p{t), 6{t), 'y{t). It can be determined from second equation (4.33) as 
a result of simple integration: 



(4.42) 7(i) = j 



* vit) sinieit)) , 

^ ' ^ ^ " dt + const . 



Formulas (4.37), (4.39), (4.41), and (4.42) prove that Newtonian dynamical system 
with force field (4.31) is integrable in quadratures. 

§ 5. Special classes of solutions of normality equations. 

Non-trivial examples of dynamical systems admitting the normal shift, which 
are considered in §4, were found in paper [Bol3] (see also preprint [Bol2]). Their 
existence stimulated further investigations and caused the rise of the theory of New- 
tonian dynamical systems admitting the normal shift on arbitrary Riemannian and 
Finslerian manifolds. Thereby some special classes of such systems were found. 

1. Geodesic flows of conformally euclidean metrics. Transferring classi- 
cal Bonnet construction to Riemannian geometry one obtains the construction of 
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geodesic normal shift, which reproduce all properties of initial construction. In par- 
ticular this means that geodesic flow is a Newtonian dynamical system admitting 
the normal shift of hypersurfaces in the metric of that manifold, where it is deter- 
mined. The space M = turns to Riemannian manifold, if we define conformally 
Euclidean metric in it with the following components: 



(5.1) Qij = e '^^Sij, where 4 



1 for i = j, 
for i ^ j. 



Here / = /(r) is some scalar field determining conformal factor e~'^^ in metric (5.1). 
Geodesic flow of conformally Euclidean metric (5.1) in standard metric of is 
represented by Newtonian dynamical system with force field 

(5.2) F = -|v|2.V/ + 2(V/, v)-v. 

Conformal factor e~'^^ changes lengths of vectors when transferring from standard 
Euclidean metric in to the metric (5.1). However, it doesn't change the measure 
of angles between vectors. For this reason force fields (5.2) admit the normal shift 
of curves in standard Euclidean metric in M?. 

2. Metrizable dynamical systems. Problem of metrizability arises in con- 
nection with force fields of the form (5.2). The matter is that the set of trajectories 
of one Newtonian dynamical system with force field Fi can be immersed into the 
set of trajectories of other dynamical system with the force field F2. In such situa- 
tion we say that first system inherits trajectories of second one. Dynamical systems 
that inherit trajectories of the system with force field (5.2) were called metrizable 
systems. The problem of metrizability consists in describing all m,etrizahle Newto- 
nian dynamical systems admitting the normal shift. This problem was investigated 
and was completely solved in papers [Shrl] and [Shr3]. Force fields of metrizable 
dynamical systems admitting the normal shift of curves in M? have the following 
form 

(5.3) F(r, v) = -|v|2 • V/ + 2 (V/, v) • v + ^ • H{\w\ e'^ e^. 

Force fields (5.3) are parameterized by two smooth scalar functions: the function of 
two variables / = f{x,y) and the function of one variable H = H{v). 

3. Dynamical systems of multidimensional type. Force fields of the form 
(5.3) are trivial in some sense: normal shift of curves along trajectories of Newtonian 
dynamical systems with such force fields coincides with geodesic normal shift defined 
by metric of the form (5.1). Therefore the problem of constructing non- metrizable 
dynamical systems admitting the normal shift was actual. In two-dimensional n = 2 
case this problem is solved by the examples 1 and 2 considered in § 4. Some examples 
of non-metrizable dynamical systems admitting the normal shift in multidimensional 
case n ^ 3 were constructed in paper [Brol]. Then we had the problem of exhaustive 
description of all solutions of complete system of normality equations in multidimen- 
sional case. This problem was considered in paper [Bol9]. However, there an error 
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in calculations were committed. This error was corrected in [Shr5] . As a result the 
following theorem was proved. 

Theorem 5.1. Newtonian dynamical system in Riemannian manifold M of the 
dimension n ^ 3 admits normal shift if and if its force field F is given by formula 

where W is arbitrary function on M depending on additional numeric parameter 
V = |v| and satisfying the condition Wy = dW/dv ^ 0, while h = h{v) is an 
arbitrary function of one variable. 

In two-dimensional case n = 2 force fields of the form (5.4) also satisfy the normal- 
ity equations (see theorem 5.2 below). Corresponding dynamical systems in M = 
are called systems of multidimensional type. 

Theorem 5.2. Force fields of the form (5.4) satisfy the normality equations 
(1.16) and (1.17). 

Proof. The equations (1.16) and (1.17) are equivalent to the equations (1.1), 
which they were derived from. Therefore we can consider the expansion (1.3) for F 
and check the equations (1.1) for A and B. Coefficients A and B in the expansion 
(1.3) are determined by formulas (1.5) and (1.7). For the field (5.4) they yield 

The function W ~ W{r, v) depends on velocity vector v only through its dependence 
on = |v|. Therefore \7W = Wy ■ N. Using this fact, we can explicitly calculate the 
velocity gradient for the function A: 

Then we find the coefficients as and a4 in the expansions (1.4): 

(VW, M) 



(5.6) ai 



Wy 



(VW^, v) h{W) (VWy, v) (VW, N) 



Let's compare (5.6) with formula for B in (5.5) and get sure that first normality 
equation (1.1) for the force field (5.4) is fulfilled. 
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Now let's check second normality equation in (1.1). First let's calculate velocity 
gradient of the function B: 

VB = . N + iZl^ . . _ , . _ (Z^ . M. 

Wy Wy {W^y Wy 

Then we find coefficients (3^ and fii in the expansions (1.4): 
(5.8) Pa = - ' 



Wy 



Coefficients a2 and /3i in the expansions (1.4) are determined as follows: 

a2 = (VA M>, A = {WB, N). 

By direct calculations according to these formulas we get: 

(5.10) 

VmVvW^ ^ (VW, v) {VWy, M) 



(5.11) /3i 



Wy ■ {Wyf 

VvVMir (vir. M) (vir,,, v) 



Wy {Wy)^ 

Here for the sake of brevity we have introduced the following notations: 

2 2 

Vm VvVF = VvVm W = ^ ^ u> ViVjW. 

i=i j=i 

Adding the equalities (5.10) and (5.11) for the sum a2 + (3\ we obtain formula 



- + ''-^(™'.M>-;ff|(V«'.,M> 
(5.12) 

{WW, v) {WWy, M) _ (VW, M) {WWy, v) 

^ (wg^ Wv? ■ 

Now let's calculate as — /34. Taking into account (5.7) and (5.8), we find 

NW,v)-hiW) ^ (VW„,v) 

(5.13) as - /34 = ^ + (W^) - ^ Wy ■ 
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Using formula (5.9), we calculate another expression 
(5.14) ^ - /33 = |v| W.^ - |v| 

Let's multiply the expression (5.14) by |v|~^A, then multiply (5.13) by |v|~^B, 
and afterwards, let's add the obtained expressions and subtract (5.12) from the 
sum. Thereby let's use formulas (5.5). Left hand side of resulting equality coincides 
with left hand side of second normality equation in (1.1). While right hand side is 
identically zero. Theorem 5.2 is proved. □ 

Class of force fields of multidimensional type (5.4) comprises all previous exam- 
ples of force fields of two-dimensional Newtonian dynamical systems admitting the 
normal shift. Let's do the following substitution into the formula (5.4): 



(5.15) W{y, v)=v e-f^^K h{v) = H{v). 

As a result of substituting (5.15) into (5.15) this formula takes the form (5.3). Fur- 
ther choice H(y) = reduces it to the form (5.2). 

Examples 1 and 2 considered in § 4, are also obtained by reduction from formula 
(5.4). Let's determine a{v) as the solution of ordinary differential equation 

(5.16) a{v) = where A{v) ^ 0. 

Then we define the function /(r) = (m, r), where m = const, and we do the 
following substitution into the formula (5.4) for F: 



(5.17) W(r, v) = a{v) e'^^, h{v) = 0. 

The substitution (5.17) transforms (5.4) to the form (4.8). This corresponds to the 
example 1 in § 4 above. 

In order to obtain the force field (4.13) from the example 2 in § 4 we should choose 
another function / = /(r) in (5.17). Let's take /(r) = |r|. The function a{v) will be 
the solution of diff'erential equation (5.16) as before. 

Functions of the form (5.4) exhaust force fields of dynamical systems admitting 
the normal shift in arbitrary Riemannian manifolds of the dimension n ^ 3. Is it 
true in two-dimensional case n = 2 ? This should be clear from further analysis of 
normality equations below. 

§ 6. Point symmetries of reduced normality equations. 

Symmetry analysis is one of the powerful tools for constructing special solutions of 
differential equations. Let's apply it in order to construct new classes of Newtonian 
dynamical systems admitting the normal shift in IR.^. By means of scalar ansatz (2.4) 
we reduced normality equations (1.16) and (1.17) to one partial differential equation. 
This is reduced normality equation, which can be written as (3.18), (4.5), or (4.25). 
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All these forms of reduced normality equation arc equivalent to each other. For the 
further analysis we choose the equation (4.5) as the most simple: 

{Ay - Ae^) cos 6' - {A^ + Aey) sin 6* + 

AAe Ae Ape Ap Ay ^ AAp^ 

V V 

Let's consider the space = M** © M, coordinates of points in which are the inde- 
pendent variables x, y, v, 9 of the equation (6.1) and the variable A. Here A is also 
treated as independent variable: 



(6.1) 



(6.2) 



e TMi 



Vector field U in five-dimensional space from (6.2) can be defined as linear dif- 
ferential operator of the first order 

d d d d d 



Coefficients X, Y, V, Q, and a in differential operator (6.3) are the functions of x, 
y, V, 6, A. Vector field (6.3) corresponds to the local one-parametric group of local 
diffeomorphisms in (see [Kobl]): 



X 




X 




x{x,y,v,9,A,t) 


y 




y 




y{x,y,v,e,A,t) 


V 


1 — > 


V 




v{x,y,v,e,A,t) 












0{x,y,v,9,A,t) 


A 




A 




A{x,y,v,0,A,t) 



Functions x, y, v, 0, A determining the transformation (6.4) should satisfy the 
system of differential equations with respect to parameter t: 



f ~ 
X = 


X{x 


,y 


V 


,e,A) 


y = 


Y{x 


y, 


V, 


e,A), 


V 


V{x 


y, 


V, 


9, A), 


II 


Q{x, 


y, 


V, 


e,A), 


/a = 


a{x, 


y, 


V, 


9, A). 
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For t = the map (ft in (6.4) is the identical map. This expressed by the conditions 
(6.6) 



L=o 



y\ 



lt=0 



y, 



lt=0 



lt=0 



A\ = A. 



Assuming vector field U in (6.3) to be smooth, we consider Taylor expansions of the 

functions x, y, v, 0, A from (6.4) with respect to parameter t in the neighborhood of 
the point t = 0. Thereby we restrict ourselves by the terms of the first order. From 
the relationships (6.5) and (6.6) we derive: 



(6.7) 



x = x + X{x, y, v,0,A)-t + 0{t^), 
y = y + Y{x,y,v,e,A)-t + 0{t^), 
V = v + Vix, y, V, e,A)-t + 0{t^), 
e ^ d + e{x,y,v,9,A) -t + Oit^), 
A = A + a{x, y, v,9,A)-t + 0{f). 



In order to invert the transformations forming one-parametric group in t we should 



only change t by —t, i e. {(pt) ^ 
can be written as follows: 



— (fi-t- Therefore transformations inverse to (6.7) 



(6.8) 



X 




X 




x{x,y,v,e,A,-t) 


V 




y 




y{x,y,v,9,A,-t) 


V 


1 — > 


V 




v{x,y,v,9,A,-t) 


9 




9 




9{x,y,v,9,A,-t) 


A 




A 




A{x,y,v,9,A,-t) 



The expansions analogous to (6.7) for transformations (6.8) are written as 



(6.9) 



x = x- X{x, y, v,e,A)-t + 0{t^), 
y = y-Y{x,y,v,e,A)-t + 0{t^), 
v = v- V{x, y, v,e,A)-t + 0{t^), 
6 = 6- Q{x, y, v,6,A)-t + 0{t^), 
A = A- a{x, y, v,6,A)-t + 0{t'^). 



Each function A = A{x, y, v, 6) determines some hypersurface F in five-dimensio- 
nal space (6.2), F is a graph of the fimction A = A{x,y,v,9). The solution of the 
equation (6.1) is exactly that function we need. Applying transformation (6.7), we 
get new hypersurface Ft = </'i(F). For sufficiently small value of parameter t — > 
hypersurface Fj slightly differs from initial hypersurface F. The hypersurface Fj or 
some part of this hypersurface can be considered as a graph for some other function 
A = A*{x,y,v,6). So the transformations (6.4) defined by vector field (6.3) can be 
treated as transformations in the set of functions A = A{x, y, v, 6). 
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Definition 6.1. Vector field (6.3) is called the field of point symmetry for the 
equation (6.1) if corresponding transformations ipt in (6.4) transforms each solution 
of this equation into another solution of the same equation. 

Field of point symmetry of the equation (6.1) form Lie subalgebra in the Lie algebra 
of all vector fields in five-dimensional space (6.2). Denote it by £. Algebra £ is 
called the algebra of point symmetries of the equation (6.1). Presently there is a 
well-developed theory for calculating point symmetry algebra of arbitrary differential 
equation (see more details in [Ibrl] and [Olvl]). Let's apply it to the equation (6.1). 

Let A{x, y, v, 6) be the solution of the equation (6.1). Wc substitute the function 
A = A{x, y, V, 9) into the arguments of the functions x, y, v, 6, A, which determine 
the transformation (6.4). As a result we get five functions 

X = x{x, y, V, e, A{x, y, v, 0).t), 
y = y{x, y, V, e, A{x, y, V, 6), t), 
(6.10) V = v{x, y, V, e, A{x, y, v, 0),t), 

= 0{x, y, V, 0, A{x, y, v, 0),t), 
A = A{x, y, V, 0, A{x, y, v, 0), t) 



of four variables (the variable t here is considered as parameter). Functions (6.10) 
in implicit (parametric) form determine the function A = A*{x,y,v,0). In order to 
calculate this function explicitly we should use first four relationships (6.10) for to 
express x, y, v, and through x, y, v, and 9. Then the obtained expressions for 
x, y, V, and 9 should be substituted into fifth relationship (6.10). Denote by J the 
Jacoby matrix of the change of variables given by first four relationships (6.10): 

(6 11) d{x,y,v,9) 
^ ' d{x,y,v,0y 

Let's consider some components in Jacoby matrix (6.11): 

Ji = Xx{x,y,v,0, A,t) + XA{x,y,v,0, A,t) ■ A^, 
Jl = Xy{x, y, V, 0, A, t) + xa{x, y, v, 0, A, t) ■ Ay, 

(6.12) 

J3 = ey{x,y,v,9,A,t) + eA{x,y,v,0,A,t) ■ Ay, 
Jf = ee{x, y, V, 9, A, t) + 9 a [x, y, v, 9, A, t) ■ Aq. 



From (6.12) we see that components of Jacoby matrix (6.11) depend on variables x, 
y, V, 9, and A. Moreover, they depend linearly on derivatives A^, Ay, Ay, and A0: 



(6.13) 



= J^{x,y,v,9,A,A^,Ay,Ay,A0,t). 
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Denote by / = J ^ the matrix inverse to J. Its components depend on the same 
quantities as components of J in (6.13): 



(6.14) 



= I^{x, y, V, e, A, A^, Ay, A^, Ae, t). 



Though the dependence on derivatives A^, Ay, A^, and Aq in (6.14) is not hnear. 

We use components of matrix / for to calculate partial derivatives of the function 
A*{x, y, V, 9), which is defined implicitly by the relationships (6.10). Remember that 
the variables x, y, v, are expressed through x, y, v, 9 due to first four relationships 
(6.10), and / is Jacoby matrix for such expressions: 



(6.15) 



d{x,y,v,9) 
d{x,y,v,9) 



Then, due to (6.15), for the derivative A^. from (6.10) we obtain the expression 



(6.16) 



Al = {A^{x,y,v,9,A,t) + AA{x,y,v,9,A,t)-A^) -11 + 
+ {Ayix,y,v,9,A,t)+AA{x,y,v,9,A,t)-Ay) ■ + 
+ {A^{x,y,v,9,A,t) + AA{x,y,v,9,A,t)-Ay) ■ if + 
+ {Ae{x,y,v,9,A,t) + AA{x,y,v,9,A,t) ■ Ae) ■ If. 



We have similar expressions for other derivatives: 



(6.17) 



(6.18) 



(6.19) 





= {Ax{x,y,v, 


9,A,t) 


+ AA{x,y,v, 


9,A,t) 


■ Ax) 






+ {Ay{x,y,v 


9,A,t) 


+ AA{x,y,v, 


9,A,t) 


' ^v) 






+ {Av{x,y,v, 


9,A,t) 


+ AA{x,y,v, 


9,A,t) 


' A^) 


■il + 




+ {Aeix,y,v, 


9,A,t) 


+ AA{x,y,v, 


9,A,t) 


■Ae) 


■il 




= {Axix,y,v, 


9,A,t) 


+ AA{x,y,v, 


9,A,t) 


■ Ax) 






+ {Ay{x,y,v 


9,A,t) 


+ AA{x,y,v, 


9,A,t) 


■ Ay) 






+ {Ay{x,y,v 


9,A,t) 


+ AA{x,y,v, 


9,A,t) 


' Ay) 






+ iA0{x,y,v, 


9, A, t) 


+ AA{x,y,v, 


9,A,t) 


■Ae) 




Ai 


= {A^{x,y,v, 


9,A,t) 


+ AA{x,y,v, 


e, A, t) 


■ Ax) 






+ {Ay{x,y,v, 


9, A, t) 


+ AA{x,y,v, 


9, A, t) 


■ -^y) 






+ {Av{x,y,v, 


9,A,t) 


+ AA{x,y,v, 


9, A, t) 


' Ay J 


■il + 




+ {Ae{x,y,v, 


9,A,t) 


+ AA{x,y,v, 


e, A, t) 


■Ae) 


■Il 



Let's distract ourselves for a while from particular form of right hand sides of the 
equalities (6.16), (6.17), (6.18), and (6.19). Let's rewrite them in symbolic form 
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keeping only the indication of functional dependence of those quantities, which are 
present in right hand sides of these equalities: 

Ai = Ai{x,y,v,e,A,A^,Ay,Ay,A0,t), 
4 = ^2(a;, y, V, e, A, A^,Ay,Ay,A0,t), 
A% = Ai{x,y,v,e,A,A^,Ay,A^,A0,t), 
A\ = Ai{x,y,v,e,A,A:„Ay,A^,Ae,t). 

The functions Ai, A2, A3, and A4 (6.20) do not depend on particular choice of 
function A = A{x,y,v,6). The function A (its explicit form) is of importance 
only when we substitute its derivatives A^, Ay, A^, Ag into the arguments of the 
functions (6.20). If we consider A, A^, Ay, Ay, Ag as independent variables and if 
we complete x, y, v, and A in (6.4) by functions (6.20), then we obtain the extension 
of transformations (fit from to M^. It is called the first extension of local one- 
paramctric group of local transformations ipf. 

Differentiating the relationships (6.20), we obtain analogous formulas for second 
derivatives A^,-, A^-, A^,-, Ai-, A--, A*~~, A* -, A\-, A* -, A*~~: 

^Ix = ^n{x,y,v,9,A,A^,Ay,Ay,A0,A^x,--- ,Agg,t), 
A%^ = Ai2{x,y,v,6,A,A^,Ay,Ay,Ae,A^^,... ,Agg,t), 
^xv = ^13(2;, y, V, e, A, A^, Ay, Ay,Ag, A^^, ... , Agg,t), 
Aig = Ai4{x, y, V, 6, A, A^, Ay, Ay, Ag, A^^, ... , Agg,t), 

(6.21) 

Aig = A24{x, y, V, 6, A, A^, Ay, Ay, Ag, A^^, ... , Agg,t), 

Aly = As3{x,y,v,9,A,A^,Ay,Ay,Ag,A^^,... ,Agg,t), 
Al§ = As4{x,y,v,e,A,A^,Ay,Ay,Ag, A^^ ,Agg,t), 
A\g = A4i{x, y, V, 6, A, A^, Ay, Ay, Ag, A^^, ... , Agg,t). 

Ten functions (6.21) define second extension of local one-parametric group of 

transformations ipt. Second extension acts in the space M.^^. 

First and second extensions of transformations (ft form local one-parametric 
groups of local transformations in the spaces and respectively. They corre- 
spond to some vector fields in and R-*^^, these vector fields are called first and 
second extensions of vector field U from (6.3). First extension can be represented 
by the following differential operator: 

U X Y ^ V O ^ ^ 
dx dy dv 86 OA 

(6.22) 
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For the second extension we have similar differential operator, but with more terms: 

U X ^ Y 

dx dy dv 89 dA 



d d d d 



■ ai —. 1- 02 TT-. 1- CIS — h 04 



(6.23) 



dAy dA,, dAg 
+ ^ + ^^4- ^ + ^ 

dAxx dAxy dAxv dAx6 

d d d d d 

+ 022 TT-A 1" <*23 TT-A ^ ^^24 TT-a ^ <l33 TT-a ^ "^34 771 1" 044 



dAyy dAyv 8 Ayg OAyy OAyg 8 AqO 

Coefficients ai, a^, 03, 04, an, 012, 013, 014, 022, 023, «24, 033! 034, and 044 in (6.22) 
and in (6.23) can be calculated in explicit form. For to calculate a\, 02, 03, 04 one 
should consider Taylor expansions of the functions (6.17), (6.18), (6.19), and (6.20) 
with respect to parameter t in the neighborhood of the point t = Q: 



A% = Ax + ai{x,y,v,e,A,Ax,Ay,A^,Ae)-t + 0{t^), 
Al = Ay + a^ix, y, v, 0, A, Ax, Ay, Ay,Ae)-t + 0{t^), 

A% = A.,+ai{x,y,v,e,A,Ax,Ay,A^,Ag) -t + Oit^), 
A\ = Ax + ai{x,y,v,e,A,Ax,Ay,A^,A0)-t + O{t^). 



Detailed calculations leading to (6.24) are standard. They are not of particular 
interest here. We shall give only the result of these calculations. For the function 

ai{x,y,v,9,A,Ax,Ay,Ay,Ag) we get 

ai=ax+ aA Ax - Ax Xx - Ay Yx - Ay Vx - Ag 9^- 

(6.25) 

-AxXaAx- AyYAAx- A^VaAx- AsQaAx. 
Similar expressions are available for 02, 03, and 04: 

(6.26) 



02 = fly + Ay - AxXy - Ay Yy - A,. Vy - Ag Qy- 

-Ax Xa Ay - Ay Ya Ay ^ Ay Va Ay - Ag Qa Ay. 



(6.27) 

(6.28) 



03 = + a A A,, - Ax X, - Ay Y„ - A„ V., - 6^- 
-Ax Xa A^ - Ay Ya A„ - Ay Va Ay - Ae Qa Ay. 

03 = + CLA Ag - AxXg- Ay Ye - AyVg - Ag Qg- 
-Ax Xa Ag - Ay Ya Ag - Ay Va Ag - Ag Ga Ag. 



Formulas for coefficients an, 0,12, ais, au, 0,22, 0,2s, 024, 033, 034, and 044 in (6.23) 
are more huge. In order to get them one should consider Taylor expansion of the 



§6. POINT SYMMETRIES OF REDUCED NORMALITY EQUATIONS. 



67 



functions (6.21) at the point t = and should take Hnear terms in such expansions: 

= ^xx + aii{x, y,v,0,A,... , A^e,Ag0) ■ t + 0{f), 
Al^ = A^y + ai2{x, y,v,0,A,... , A-,e,Aee) ■ t + 0{t^), 

(6.29) 

Ai^ = A^e + aii{x, y,v,e,A,... , A^g, Aee) ■ t + 0{t^), 
Alg = Aee + a34{x, y,v,9,A,... , A^e, Agg) ■ t + 0{t^). 

We do not give exphcit formulas for aij{x, y,v,0, A, . . . , A^g, Agg)^ they are rather 
huge. We give only the complete list of arguments in them: 

2/> ^) ^; A, Ax, Ay, A^, Ag, Axx, Axy, 

(6.30) 

AxV^ Ax0, Ayy, Ayy, Ayg , Ayy , A^g , Agg. 

Let A{x,y,v,9) bo an arbitrary scalar function of four arguments. Fixing some 
values of its arguments, we calculate its value and the values of all its derivatives 
listed in (6.30) for these fixed arguments. As a result we get the point of the 
space M.^^. Each point of the space R-'^^ can be obtained in this way. However, if 
choosing the function A{x, y, v, 9) we restrict our choice to the set of solutions of the 
normality equation, we would obtain only the points of some hypersurface in M}^ . 
The equation (6.1) itself is the equation of such hypersurface provided that we treat 
all quantities in (6.30) as independent variables. Remember that reduced normality 
equation equation (6.1) has the following structure: 

(6.31) F{x,y,... ,Ayg,Aee) = Q. 

Function F in (6.31) is determined by left hand side of the equation (6.1). According 
to general theory of point symmetries (see [Ibrl] and [Olvl]) vector field (6.3) is a 
field of point symmetry for differential equation (6.31), if its second extension is 
tangent to hypersurface in M^^ defined by the equation (6.31). This condition is 
written as 

^,dF ,,dF ^BF ^dF OF OF OF 
' = ''d^^'^d^^''^^^m^^dA+^^dAx^^^dAy^ 

OF dF dF dF dF OF 

(6.32) +03 iT-r + 04 ^T-r + an — — + jr-r— + ai3 tt-t— + au 



dAy dAg dAxx dAxy dA^v dA^g 

OF dF OF dF dF dF 

+ 0,11 -TT-, 1" ^123 771 1- a24 -TT^ ^ <l33 -7r\ 1" ^34 TTa ^ <l44 



dAyy dAyy d Ayg dAyy BAyg dAgg 

The equality (6.32) should be fulfilled at all points of the space M^^, where the 
equality (6.31) is fulfilled. The equality (6.32) is called the determining equation 
for the field of point symmetry. 
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§ 7. Calculation of the point symmetries. 

Let's consider the determining equation (6.32) for the field of point symmetry U 
applied to our case, when the equation (6.31) has the form (6.1). The equation (6.1) 
can be resolved with respect to one of the second order derivatives: 

V {Ay - Aq:^) cos 9 AgAy 
Aye = J + ^ 

(7.1) 

_ v{A^ + Aey) sine _^ ^ _j. Age 
A V vA 

Substituting (7.1) into (6.32), we obtain the determining equation for the field of 
point symmetry U in the following form: 

(7.2) Fdet{x, y, . . . , Ayy,Agg) = 0. 

In left hand side of (7.2) we have all variables (6.30) except for Ayg The equality 
(6.32) brought to the form (7.2) should be fulfilled identically in the space R^^. 

The equation (7.1) is polynomial in derivatives. Therefore the variables A^, Ay, 
Ay, Ag, Axxi Axy, Axy, Aj;0 , Ayy , Ayy, Ayg, Ayy , Agg euter the equality (7.2) 
polynomially. The entries of variables A^x, Ayy, Ayy, and Agg in the equality (7.2) 
are linear. In highest order with respect to A^x the equation (7.2) looks like 

(7.3) A C0s9{Xg +XAAg)-Axx + -.- = 0. 

Here Ag is independent variable. The quantities v, A, cos in this equation cannot be 
zero simultaneously. Therefore from the equality (7.3) we obtain that the following 
two partial derivatives are zero: 

(7.4, X. = f=0, ^^ = §J = ». 

This means that first component X in the field of point symmetry (6.3) shouldn't 
depend on the variables 9 and A. Analyzing the entries of derivatives Ayy and Ayy 
in (7.2) we obtain analogous results: 

dY dY 

(7.5) ^^ = ^=0. ^-=^4=°' 

dV dV 

(7.6) ye = ^ = 0, y^-0A='- 

Taking into account the relationships (7.4), (7.5), and (7.6), wc can split the equation 
(7.2) with respect to the variable Agg, the entry of which is hnear: 

^det(a=, • . • , Ayg) ■ Agg + F^l{x, Ayg) = 0. 
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Prom this equation we obtain two separate equations. We write them as follows: 

-^det(^'2^'^' • ■ ■ ,Ayv,Ay0) = 0, 

(7.7) 

Calculating _FjJ,[ explicitly, wc find that first equation (7.7) doesn't contain actual 
dependence on A^y, A^v, A^e, Ayy, Ayg. While with respect to Ag left hand side of 
this equation is a polynomial of the second order: 

(7.8) -2AveA- Ae + ... = 0. 

Prom the relationship (7.8) we see that the partial derivative 9^ should vanish: 

(99 

(7.9) ^-=0A='- 

This means that fourth component of the vector field of point symmetry (6.3) do 
not depend on A. With the equality (7.9) in mind, the equation (7.8) turns to be 
linear in Ag. We write it as follows: 

F:^:,\x,...,A,).Ae + F!^:,\x,...,A,) = 0. 

The above equation splits into two separate equations: 

(7.10) i^r=o, i^r=o. 

The equations (7.10) are sufficiently small. We can write them expficitly. These 
equations are written as follows: 

K + Vy sin^e + ^ = 

(7.11) 

AV a AttA ^AQe 
= ^ + ^+2 -. 

V V V V 

AO 

(7.12) 9^ cos^e + 9„ siri^e + = -ae- 

V 

Now let's consider the second equation (7.7). It is linear with respect to variables 
Axd, Ayg. We write it in the following form: 

(7.13) F^°:t' ■ + FZ' ■ Aye + = 0. 

The equation (7.13) splits into three separate equations: fjet^ = 0, i^jgf = 0, and 
-^det' = 0. The last equation doesn't contain second order derivatives 

-^xyf -^xvi -^xBt 
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Ayy, Ay0. Lcft liaiid side of this equation is a polynomial of the third order with 
respect to variable A0 . In the leading term we have 

(7.14) vAaAA-A0+... = O. 

Prom the relationship (7.14) we get gaa = 0. This means that fifth component of 
the vector field of point symmetry (6.3) is linear in A: 

(7.15) a = a{x, y, v, 6) + /3(.t, y, v, 9) ■ A. 
Moreover, from the relationships (7.4), (7.5), (7.6), and (7.9) we get 

X = X{x,y,v), Y = Y{x,y,v), 

(7.17) 

V = V{x,y,v), e = eix,y,v,9). 

Thus, the entries of the variable A in components of the field of point symmetry 
(6.3) are calculated in explicit form. 

Let's substitute (7.15) into the equation (7.11). As a result we obtain the equation, 
which is linear in A. It splits into two separate equations: 

V 

(7.18) Vy = 29^, a = cos6' + wK sin6'. 

V 

Second equation (7.18) expresses the fimction a{x,y,v,6) through Vx and Vy in ex- 
plicit form. Let's substitute this expression into (7.15), then substitute the resulting 
expression for a into the equations F^J^'' = and F^J^" =0. As a result of such 
substitution left hand sides of these equations turn to be polynomials of the second 
order with respect to A. In leading order they have the following form: 

(7.19) XyA^ + ... = 0, YyA'^ + ... = 0. 

From (7.19) we get X„ = and Yy = 0. Now we can specify the relationships (7.17): 



(7.20) X = X{x,y), Y = Y{x,y), 

Due to (7.20) in the equations F^^^ = and F^^^ = we can keep only terms 
linear in A, zero order terms in A being identically zero. First order terms in these 
equations can be written explicitly. Thus we can reduce (7.19) to the equations that 
do not contain the variable A at all: 

{Xy + @) sin6» 

(7.21) 

{Y^-Q) cos 6* 



Vy + ^ - P-Yy]sme. 
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Let's multiply the equation (7.21) by sinO, multiply second one by cos^^, then let's 
subtract second equation from the first one. At a result we obtain the equation 
expressing parameter O through X^, Xy, Y^, and Yy-. 

(7.22) e = cos^e - Xy sm^e + {Yy - X^) sin 6* cos 6*. 

Further let's consider again the equations (7.11) and (7.12). First of them is 
already reduced to the pair of separate equations (7.18), one of which expresses 
a through T4 and Vy. Let's substitute this expression for a and the expression 

(7.22) for 8 into the first equation (7.18). Upon rather simple transformations of 
trigonometric functions we get the equation 

(7.23) [Yy - X^) cos{2e) - [Y^ + Xy) sin(2^) = y - ^- 

The variable 9 enter to the equation (7.23) only through trigonometric functions 
cos(20) and sin(20) (see the relationships (7.17) and (7.20) above). Therefore (7.23) 
splits into three separate equations. One of them is the following 

(7.24) K = -• 

V 

Other two equations (it's remarkable) have the form of Cauchy-Riemann equations 
with respect to functions X{x,y) and Y{x,y): 

(7.25) Xy = —Yx, Xx = Yy. 

Let's compare (7.24) and (7.18). From such comparison we get Qg = 0. The 
equation (7.24) itself is, in essential, an ordinary differential equation for V. Its 
general solution is given by V = w W{x, y). Now we can specify (7.17), writing these 
relationships in the following form: 

X = X{x,y), Y = Y{x,y), 

(7.26) 

V = vW{x,y), Q = e{x,y,v). 

Note that the last relationship (7.26) doesn't contradict to the formula (7.22) for 9. 
Indeed, if we take into account Cauchy-Riemann equations (7.25) and use trigono- 
metric identity s\v?9 = 1 — cos^0, then we get 

(7.27) Q = -Xy= Y^. 

Formula (7.27) gives further specification of the relationships (7.26): 

X = X{x,y), Y = Y{x,y), 

(7.28) 

V = vW{x,y), Q = e{x,y). 
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In the next stop wc substitute (7.27) into any one of the relationships (7.21). The 
results of these two substitutions coincide with each other. They yield the equality 
determining parameter /? in (7.15): 

(7.29) l3 = 2W-Yy. 

Prom the whole set of the above equations now remain two equations. These arc 
the equation (7.12) and the equation = 0, which was derived from (7.13). 

Let's substitute the reduced functions (7.28) into the equation (7.12). Thereby we 
take into account formulas (7.18) and (7.29) for a and /?, formula (7.27) for Q, and 
formula (7.15) for parameter a. We keep in mind Cauchy-Riemann equations (7.25) 
for X and Y as well. Then we obtain the equation 

(7.30) {Wy - X^y) sin e + {X^^ - W^) cos ^ = 0. 

The equation (7.30) splits into two separate equations, which form the system of 
Pfaff equations with respect to the function W = W{x,y): 



(7.31) 




System of Pfaff equations (7.31) is compatible due to the Cauchy-Riemann equations 
(7.25). It's easy to integrate it. As a result we determine the function W{x,y) up 
to some arbitrary constant C: 

(7.32) W = X^ + C. 

Now it remains to consider only one equation F^°^ =0. In a highest order in A 
it yields uaa = in (7.14), which, in turn, yields (7.15). In lower orders in A this 
equation appears to be identically fulfilled due to above formulas for the functions 
X, y, V, e, and o. 

§ 8. Algebra of point symmetries of the reduced normality equation. 

Results obtained in § 7 allow to give complete description of the structure of any 
vector field of point symmetry for the equation (6.1). It's convenient to do it in 
terms of complex variable z — x + iy. Due to Cauchy-Riemann equations (7.25) 
the components X and Y in vector field (6.3) are real and imaginary parts of some 
holomorphic function ^(z). We write them as 

(8.1) X = Re^{z), y = Im^(^). 

The relationships (7.27) and (7.32) now are written as follows: 



(8.2) 



W = Re^'{z) + C. 
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In addition to z ^ x + iy wc introduce another complex variable u = u e'^ . Then 
from the relationships (7.18) and (7.29) we get 

(8.3) a = Re{uf{z)), (3 = Re^'{z) + 2C. 

Let's substitute (8.3) into the formula (7.15). Thereafter we can write formula for 
vector field of point symmetry 



2 dA dA 

Theorem 8.1. Each vector field U of point symmeti-y of reduced normality equa- 
tion (6.1) is determined by some holomorphic function ^{z) and some real constant 
C according to the formula (8.4). 

Let U\ and U2 be two vector fields of point symmetry for the equation (6.1). 
Suppose that first of them is determined by the function ^i(z) and by constant C\, 
while second is determined by function ^2{z) and by constant C-j- By means of direct 
calculations we find that commutator of the fields C/3 = \U\, U^] has the form (8.4). 
It is determined by holomorphic function 

(8.5) U^) = ii{z)^^{z)-U^)^[{z) 

and by constant C3, where C3 = 0. We can associate functions ^i(^) and ^2 (-2) with 
the following two holomorphic vector fields: 

(8.6) 6(.)|, 6(.)|. 

Then function (8.5) is associated with commutator of holomorphic vector fields (8.6). 
Adding constants Ci and C2 to vector fields (8.6) corresponds to the construction 
known as central extension (see [Brbl], Chapter I, § 1, subsection 7). 

Theorem 8.2. Lie algebra £ of point symmetries of normality equation (6.1) is 
isomorphic to central extension of Lie algebra of holomorphic vector fields in<C = M?. 

§ 9. Invariant solutions. 

Let U be some vector field of point symmetry of reduced normality equation 
(6.1). It generates some one-parametric group of point transformations (6.4). Say 
that function A{x,y,v,6) is invariant with respect to vector field U in (6.3) if 
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function A''^{x,y,v,9) defined implicitly by the relationships (6.f0) do not depend 
on t and coincides with A{x, y, v, 9). This condition is expressed by the relationship 



(9.1) 



A = A{x,y,v,e) 



which binds functions x, y, v, 9, and A in (6.10). Let's differentiate the relationship 
(9.1) with respect to parameter t. Then substitute t = and take into account the 
expansions (6.7). As a result we get the relationship 



This relationship expresses the condition of inv£iriance of the function A{x, y, v, 9) 
with respect to vector field (6.3) in the form that can be checked by direct calcula- 
tions. Briefly this relationship is written as 



Definition 9.1. The solution A{x,y,v,9) of reduced normality equation (6.1) 
is called invariant solution if it is invariant with respect to one or several fields of 
point symmetry of this equation. 

Vector fields of point symmetry of the equation (6.1), with respect to which some 
fixed solution A{x, y, v, 9) of this equation is invariant, form some Lie subalgcbra £,a 
in Lie algebra £ of all point symmetries of this equation. All examples of dynamical 
systems admitting the normal shift, which were considered in § 4, correspond to some 
invariant solutions of reduced normality ('filiation. Now we consider these examples 
again and calculate Lie subalgebras corresponding to them. 

Example 1. Spatially homogeneous, but not isotropic force field. Func- 
tion A{x,y,v,9) in this example has the form A = A{v) cos 9. Let's substitute 
this function into the equation (9.2). Thereby we take into account the following 
formulas for the components of the field of point symmetry: 



A^{x, y, V, 9) ■ X{x, y, v, 9, A{x, y, v, 9)) + 

+ Ay{x, y, V, 9) ■ Y{x, y, v, 9, A{x, y, v, 9)) + 
+ A^{x, y, V, 0) ■ V{x, y, v, 0, A{x, y, v, 9)) + 

+ A0{x, y, V, 0) ■ Q{x, y, v, 0, A{x, y, v, 6)) = 
= a{x,y,v,9,A{x,y,v,0)). 



(9.2) 



A^-X + Ay-Y + A^-V + As-Q = a. 



V = {X^ + C)v 



(9.3) 



e = -x, 



a = v^ {X^^ cos 9 + X^y sin 6*) + {X^ + 2 C) A. 



These formulas are obtained from the relationships (7.15), (7.18), (7.27), (7.28), 
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(7.29), and (7.32), which were derived in § 7. Substituting (9.3) and A = A{v) cos6 
into the equation (9.2), we obtain the equation: 

{X^^ cos 9 + X^y sin 6) + (X^, + 2 C) A{v) cos 6 = 

(9.4) 

= V A'{v) {X^ + C) cos e + A{v) Xy sin 0. 

By collecting coefficients of cos 61 and sin 61 in (9.4) we split this equation into two 
separate equations. They are written as follows: 

(9.5) X^yv'^ =A{v)Xy, 

(9.6) X^^ v"" + {X^+2 C) A{v) = V A'{v) {X^ + C). 

Considering the equations (9.5) and (9.6) wc find three subcases. 

Subcase 1. Function A{v) is not an exponential function, i. e. A{v) ^ k ■ v", 
where k and a are constants. In this case from (9.5) we get Xy = 0. Then from 
Cauchy-Ricmann equations (7.25) we derive = and X^ = Yy = Ci = const. 
Hence for X and Y we have the expressions 



(9.7) X = Cix + C2, Y = Ciy + C3. 

Substituting (9.7) into (9.6) and taking into accovmt that A{v) ^ c-w", we find that 
C = Ci = 0. Constants C2 and C3 remain undetermined. This means that Lie 
subalgebra £a for the solution A = A{v) cos9 in this case is two-dimensional. It is 
generated by the following two vector fields: 

= = ^y■ 

Subcase 2. Function A{v) is exponential function A{v) = k ■ v°' , where k and a 
are constants, but a ^ 2. In this case, as in previous one, the equation (9.4) splits 
into two separate equations (9.5) and (9.6). From (9.5) we derive (9.7). However, 
the result of substituting (9.7) into (9.6) is different: 

(9.9) Ci(l-a)=C(a-2). 

The equation (9.9) is solvable in form of the following two relationships: 



(9.10) Ci = (a-2)C4, C=(l-a)C4. 

Formulas (9.7) and (9.10) contain three arbitrary constants C2, C3, C4. This indicate 
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that algebra £a is three-dimensional. Here are the generators of this algebra: 



d d 
ox oy 



(9.11) 

f d d \ d d 

Subcase 3. Function A{v) is quadratic, i. e. A{v) = k ■ v"^. The equations (9.5) 
and (9.6) in this case are written as follows: 

(9.12) = kXy, = kXx- 

The equations (9.12) should be completed by the equations of Cauchy-Riemann 
(7.25). Upon rather simple analysis we can write general solution to the system of 
equations for X and Y that arises in this case: 

X = Ci cos{ky) 6*=^ + Ca sm{ky) e'=^ + C3, 

Y = C2 cosiky) e^"" - Ci sm{ky) e*^^ + C4. 

Elements of algebra £,a depend on five arbitrary constants Ci, C2, C3, C4, and C. 
Here are the generators of this algebra: 

ox oy 

d d 

1/3 = V- 2kv'^ cos6—, 

ov oA 

(9.14) U4 = e'"= cos(ky)( + kv + 2k^ v"^ cos0 + 

V ox ov oA J 



e- -n(M(| + fc|-fc^.^ sin^A 



/ d d d 

U5 = e*^"^ sm{ky) I — + kv — +2k'^ v"^ cos0 — 
\ ox ov oA 

-e'=^ cos(M(^ + fc|-fc^-^ sin^^ 

Example 2. Spatially non- homogeneous force field with marked point. 

Function A determining force field F in this case is given by formula (4.29): 

(9.15) ^^ A{v)cose ^ 
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However, this function is written in variables p, 7, v, 9. They are different from 
variables x, y, v, 9, in which the reduced normality equation (6.1) is written. One 
easily find the relations between these two sets of variables by comparing Fig. 4.1 and 
Fig. 4.2. The variable v denotes the modulus of velocity in both sets of variables. 
Angular variables differ from each other by the value of the angle 7. Therefore 
function (9.15) should be rewritten as follows: 

_ ^i^) cos(^ - 7) 



In order to complete the transfer to the variables .t, y, v, 9 one should apply the 
identity cos(6' — 7) = cos 9 cos 7 + sm9 sin 7 and one should use the relationships 
binding p and 7 with variables x and y: 



X = p cos 7, y = P sm7. 

As a result one get the following expression for A: 

X cos 9 + y sinO 



(9.16) A = A{v) 



x^ + y^ 



Let's substitute (9.16) into the equation (9.2) and let's take into account the rela- 
tionships (9.3). Instead of (9.5) and (9.6) we get two other equations 



(9.17) 



P 

-2x^-^)+vA'iv)x^=0, 
P^ J P 



(9.18) 



xX + yY\ X, + C 
- 2y J +vA {v)y = 0. 



Here p — ^ x'^ + J/^- The equations (9.17) and (9.18) determine three subcases. 

Subcase 1. Function A{y) is not an exponential function, i. e. A(v) ^ k ■ v", 
where k and a are constants. In this case each of the above two equation (9.17) 
and (9.18) splits into three separate equations. Among six resulting equations we 
have two coincident ones. So it remains five equations that follows from (9.17) and 
(9.18). Here are three of these five equations: 



(9.19) X^x = 0, Xxy = 0, Xx = -C. 



78 III. ANALYSIS OF NORMALITY EQUATIONS IN R'^. 

Other two equations are more complicated. They have the following form: 

{X --xC -yXy){x^ +y^) = 2x{xX + yY), 

(9.20) 

{Y - yC + xXy){x^ + y^) = 2y{xX + yY). 

Let's complete (9.19) with the Cauchy-Riemann equations (7.25). For X and Y, 
upon solving these equations, we get 

X = -Cx + Ciy + C2, Y = -Cy-Cix + C3. 

Substituting these expressions into (9.20) we get the equations that yield C2 = 
and C3 = 0. Constants C and Ci remain undetermined. Hence Lie subalgebra £,a 
in this case is two-dimensional. Here are its generators: 

d d ., .xcoaO + ysmO d 

(9.21) 

U -- —+X— + — 
^ ^ dx ^ dy do 

Subcase 2. Function A{v) is not an exponential function, i. e. A{v) ^ k ■ v°', 
where k and a arc constants, but a 7^ 2. From (9.17) and (9.18) in this case we 
extract four equations. Here are two of them: 

(9.22) X,, = 0, X,y = 0. 

Rest two equations are written in the following form: 

2x{xX + yY) 



(9.23) X + xa{X^+C) -x{X^ + 2C) = yXy 



x^ + y^ 
2y{xX + yY) 



(9.24) Y + ya{X^ + C)-y{X^ + 2C) = -xXy+ ^ . 

X -r y 

Let's complete the equations (9.22) with Cauchy-Riemann equations (7.25). For X 
and Y this yields the following expressions: 

X = Cix + C2y + C3, Y = Ciy-C2X + Ci. 

Substituting these expressions into (9.23) and (9.24), we obtain three equations, 
solution of which is given by the relationships 



C=-Ci, 



C3 = 0, 



C4 = 0. 
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Constants Ci and C2 remain undetermined, this means that Lie algebra £,a is two- 
dimensionah dimii^ = 2. Here are its generators: 

^ d d ]f a ^ '^'^^ 6 + y sinO d 

^ ''^ dx ^ dy ^ + y^ dA ' 

(9.25) ^ ^ 
U d ^ d ^ d 

^ ^ dx 5y do 

They have the same form as generators (9.21) in previous case. 

Subcase 3. Function A(v) is quadratic, i. e. A(v) = k - v^. The equations (9.17) 
and (9.18) in this case do not split. In order to write these equations and solve 
them we use complex form (8.4) for the vector field of point symmetry (6.3). Let's 
multiply (9.18) by i and let's subtract the result from (9.17). The resulting equation 
leads to the following differential equation for holomorphic function ^{z) in (8.4): 

(9.26) ^"(^)_^^'(^) + ^^(^) = 0. 

The space of solution of linear ordinary differential equation (9.26) over complex 
numbers is two-dimensional: 

(9.27) ^{z) = Ciz'' + C2 z. 

Its dimension over real numbers is equal to 4. Besides ^(z) formula (8.4) contain real 
constant C, the value of which is not restricted by the equation (9.26). Therefore 
the dimension of algebra £,a is equal to 5: 

d d , oxcos9 + ysin9 d 

u^^^ir + v-^-kv r— ^ — , 

ox oy x'^ + y^ oA 

U d ^ d ^ d 

^ ^ dx ^ dy do' 

U -X ^ + ^ + ^ V ^ 
^ ^ dx ^ dy 2^ dv' 

d d d d 

C/4 = Xfe — + n — + kvXk-i ^ + kYk-i ^ + 
ox oy ov oO 



+ v^k({k-l) {Xk-2 cos + Yk-2 sin 0) + k Xk-i 



X cosO + y sinO \ d 

x"' + y J a A 



U5 = -Yk — +Xk^- kvYk-i +kXk-i ^ + 
ox oy ov 06 

+ v^k ( (1 - k) {Yk-2 cosO + Xk-2 smO) ^^2j^ kYk-ij —. 
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In writing components of vector fields U4 and U5 we used the following notations: 



Xg{x,y)=Re{{x + iy)''), Yg{x, y) = lm{{x + iyf ) . 



§ 10. Symmetry reduction of reduced normality equation. 

By means of scalar ansatz (2.4) the system of normality equations (1.1) were 
reduced to one differential equation. In variable x, y, v, 6 this equation looks like 
(6.1). Further reduction of this equation in general case is impossible. However, 
there are special cases, when such reduction is possible. Below we consider some of 
these special cases. 

Example 3. Spatially homogeneous, but not isotropic force field. For 
symmetry reduction of normality equation (6.1) we need to use some subalgebra Za 
in the algebra of point symmetries £ of this equation. In this example we choose 
subalgebra 2, a with generators (9.8): 

Substituting components of vector fields (10.1) into the equation (9.2), we get Ax = Q 
and Ay = 0. This reduces normality equation (6.1) to the form (4.6): 

.^„„. AAe AeAee , Aq A, A A,,, 



V V 

We can lower the order of differential equation (10.2). Let's denote 

(10.3) b=^. 

Substituting (10.3) into (10.2), we reduce (10.2) to the following quasilinear differ- 
ential equation of the first order for the function h: 

(10.4) bbe-vh^ + b^ + b = Q. 

The equation (10.4) is solved by method of characteristics (see [Karl]). Character- 
istic lines of the equation (10.4) are parametric lines in given by the equations 

(10.5) v = -v, e = b, b = -b^- b. 

It's easy to write two functionally independent first integrals for the system of dif- 
ferential equations (10.5). These are 

(10.6) Ii=e + arctan(6), I2 = — ^ , 

where u is some constant. According to the theory of quasilinear partial differential 



§ 10. SYMMETRY REDUCTION OF REDUCED NORMALITY EQUATION. 



81 



equations of the first order (sec [Karl]), general solution of the equation (10.4) is 
given in implicit form by functional equation 

(10.7) ^(7i,/2) = 0, 

where ^ is some arbitrary function of two variables. From (10.6) and (10.7) it follows 
that the equation (10.4) is integrable in quadratures. Hence the equation (10.2) is 
also integrable in quadratures. This result was obtained in [Bol8]. 

Example 4. Spatially non- homogeneous force field with marked point. 

In this example we choose Abelian Lie algebra £,a with generators (9.21): 

^ ^ dx ^ dy dA ' ^ dx ^ dy 86' 

Substituting components of these two vector fields into the equation (9.2), we get 
the system of two differential equations: 

^ Ae + -^^-^ ^ = 0, 



(10.9) "'7' + 

+ 2 T 2 + 2 ^ 2 ^ = 0- 

-\-y^ x^ -\- y^ 

The equations (10.9) are compatible. Similar to (10.4), the equations (10.9) are 
quasilinear equations of the first order. Such equations are solved by method of 
characteristics (see [Karl]). Applying this method to (10.9), we get the theorem. 

Theorem 10.1. Solution of the system of equations (10.4) is determined by an 
arbitrary function of two variables A{v,9): 

(10.10) ^^Mv,o-i)_ 



Here p = \J x^ +2/^, while the angle 7 in formula (10.10) is determined by the 
relationships x = pcos^ and y = psin^. 

Relying upon theorem 10.1, we substitute (10.10) into the normality equation 
(6.1). This determines symmetry reduction of the equation (6.1). In this case it's 
given by the following equation for the function A = A{v,9): 

i.,r,,,. AAe AeAee , Ag A^ . AAe^ 
(10.11) — 5 1 \- {A + Agg) smO = . 

V V V 

The equation (10.11) is quite similar to (10.2). But the substitution (10.3) doesn't 
lower the order of the equation (10.11). 

Example 5. Reduction determined by algebra of conformal automor- 
phisms of a disc. In this example we consider one more way for reducing the 
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normality equation (6.1). Here we use complex form of vector field of point symme- 
try (8.4) and theorem 8.2, which binds algebra £, with algebra of holomorphic vector 
fields in C = M^. Let's consider the disc K of the radius R in complex plane: 

K = {z&C: \z\ <R}. 

Group of conformal automorphisms of the domain K has the real dimension 3. Each 
automorphism in this group is determined by one complex parameter a and one real 
parameter (p (see more details in [Shbl]): 



(10.12) 



R{z-Ra) 
R — az 



For (fi ^ and a ^ we get /(z) id(z). Let's calculate the derivatives of the 
function f{z) with respect to parameters cp, a, and a at the point (p = and a = 0: 



6(^) = 

6(^) = i 



df{z) 



df{z) 



dip 



a=0 



da 

df{z) 



, df{z) 
v=o da 

a=0 



J?2 



V3=0 
a=0 



da 



dfiz) 



^=0 * da 

a=0 



R ' 



ip=0 
a=0 



R 



Functions ^i, ^2, and ^3 determine three holomorphic vector fields, which generate 
algebra of conformal automorphism of a disk K: 



f.(=)|, 



It's easy to calculate commutators of these vector fields: 



Ki,6] = -e3, [6,^3] = 46, [e3,6] = -6- 

Let's substitute the fimctions £,i{z), ^2{z), and ^3(2) into the formula (8.4) for vector 
field of point symmetry of the equation (6.1), taking C = by each such substitution. 
This determines three fields Ui, U2, Us, which are generators of three-dimensional 
non-abelian subalgebra £a in the algebra of point symmetries of the equation (6.1): 

[Ui, U2] = -C/3, [U2, C/3] = 4 t/i, [U3, C/i] = -U2. 



In order to write explicit formulas for vector fields Ui, U2, U3 in real variables x, y, 
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V, 6 we use formulas (8.1), (8.2), (8.2), and formulas (7.15), (7.26): 



(10.14) 



(10.15) 



R dx R dy 
2vx d 2yd cosO + xA d 

^ _ 2xy d _x'^ -y^ + R^ d _^ 
^ R dx R dy 

2vy d 2x d ^v"^ sm9 + y A d 

H 1-2 — . 

R dv R 88 R OA 



Lie algebra £a with generators (10.13), (10.14), (10.15) is isomorphic to matrix 

algebra so(l,2,M). 

Vector field Ui in (10.13) coincides with the field U2 in the above example 2 (see 
formulas (9.21)). For to simplify formulas (10.13), (10.14), (10.15) we transfer from 

X, y, V, 6 to variables p, 7, v, 9 used in the example 2. In these variables vector field 
(10.13) is reduced to the operator of differentiation in 7: 

(10.16) U, = ^. 

While normality equation (6.1) turns to (4.28). The invariance condition of the 
solution of (4.28) with respect to the field (10.16) is written as A^ = 0. This means 
that A = A{p, V, 9). Let's transform vector fields (10.14) and (10.15) to the variables 
p, 7, V, 9. For vector field U2 we get 

— R? ( 8 sin 7 9 \ 2v p cos 7 8 

-|- i?^ sin 7 8 cos(0 4- 7) + COS7 8 
^ R P~dj^ R dA' 

For the field U3 we get similar expression: 

p^ — R^ { . 8 cos 7 8\ 2 tip sin 7 8 
Us = ^^-f^ — sin7 7^ -^] + 



R \ ' dp p 89 J R 8v 

p^ + R^ cos'j 8 sm{9 + ^)v'^ + Ap sinj 8 
R P~dj^ R dA' 
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Writing for A ~ A{p,v,9) the invariancc conditions with respect to U2 and U3, we 
can reduce these conditions to the form of equations 

(10.17) ^^^2pv^.ine 



(10.18) Ap 



2pvA^ - 2 7'2 cose-2pA 



" R'^-p^ 
The equation (10.17) is easily solved. This determine the dependence of A on 6: 

(10.19) A = -J^—^+A{p,v). 

Let's substitute (10.19) into (10.18). As a result we get the following equation for 
the function A{p,v) from (10.19): 

(10.20, A,-^A.,^=,. 

The equation is solved by method of characteristics (see [Karl]). General solution 
of this equation contain one arbitrary function / = f{v): 

A{p,v) = vf 



B?-p^ 



Now let's substitute the above function into (10.19). Invariant solution of the equa- 
tion (4.28) then should be found in the following form: 

. 2 pv"^ cos 9 
(10.21) A = -^-^-^+vf 



R^-p^ '\R^-p^ 

Substituting (10.21) into (4.28), wo find that (10.21) is the solution of the equation 
(4.28) for any choice of function / = f{v). In variables x, y, v, 6 for A we have 

(10.22) ^^_ 2vH.cose + y^ne) ^^^ 



_ ^2 _ y2 y]l2 _ J.2 _ y1 

Function (10.22) is the solution of normality equation (6.1). 

§ 11. System of multidimensional type 
and specifically two-dimensional systems. 

The solutions of normality equation considered in§ 9 correspond to Newtonian 
dynamical systems with force fields (4.8) and (4.31). If A{v) = k ■ v"^ these systems 
are metrizable (see formula (5.3) in §5). For other choice of A{v) they aren't 
metrizable, but they belong to more wide class of systems of multidimensional 
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type (see formula (5.4) and formulas (5.16), (5.17) in §5). Parameter A for the 
force fields (5.4) is determined by formula (5.5): 

(11.1) ^^M^)-Jw,v)_ 

In variable x, y, v, 6 formula (11.1) looks like 

_ h{W) - V {W^ cos e + Wy sin 6*) 



(11.2) 



Remember that here h = h{W) is a function of one variable, while W = W{x, y, v) 
is a function of three variables. 

In multidimensional case n ^ 3 functions (11.1) exhaust all solutions of reduced 
system of normality equations (see theorem 5.1 in §5 above and thesis [22]). Two- 
dimensional case is quite different in this point. Here there is only one reduced nor- 
mality equation. It has solutions with no multidimensional analogs. Let's con- 
struct an example of such essentially two-dimensional solution using the relationship 
(10.7) from § 10. For this purpose we choose ^(/i, 12) = (cos(/i) — 12)^ — 1/2. Then 
one can solve the equation (10.7) in explicit form and find the function b = b{v, 9): 



g ^ sin 2^ -|- 2 u u cos 9 + v^/v'^ + Auv s\n9 + 2v? 

4 u w sin 6* 4- 2 u2 _ cos 26* 

According to (10.3), function 6(w, 6*) is logarithmic derivative of the required function 
A{v,9). Hence A = e^^'"'''\ where 



(11.4) 



^ . sin 20 + 2 w u cos 9 + v 
4 w i> sin -|- 2 _ ^2 ^.Qg 20 



+ 4 ?i 7; sin + 2 
4u V sm9 + 2u'^ — v"^ cos 29 



First of two integrals (11.4) can be calculated explicitly in elementary functions. 
Second one, upon the change of variable t = sin0, is brought to elliptic integrals 
(see [Gral]) and expressed through elliptic functions (see [Ahil]). The dependence 
on theta expressed by elliptic functions cannot be obtained by formula (11.2), no 
matter what function W = W{x,y,v) is chosen in it. The variable 9 enters the 
formula (11.2) in purely trigonometric form through sin^ and cos 9. This means that 
the solution of normality equation given by integrals (11.4) is not among solutions 
of multidimensional type. It is essentially two-dimensional solution. 
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DYNAMICAL SYSTEMS ON RIEMANN SURFACES. 



§ 1. Newtonian dynamical systems in 
two-dimensional Riemannian manifolds. 

Let M be two-dimensional manifold with metric tensor g and corresponding met- 
ric connection F. Newtonian dynamical system in M is determined by a system of 
ordinary differential equations of the form 



(1.1) 



V, 



Vtv = F(x,v), 



where x is a vector composed of local coordinates x^, of a point p £ M in some 
local map, while v is a tangent vector at this point: 



The solution of the equation (1.1) determines a parametric curve p = p{t) in M, 
which describes the motion of a point of unit mass in the force field F. By VjV 
in (1.1) we have denoted the covariant derivative of velocity vector v with respect 
to parameter t along this curve. This is the vector whose components in local 
coordinates are determined by formula 



(1.2) 



i=i j=i 



1, 2. 



Trajectories of dynamical systems of the form (1.1) are completely determined by 
initial point pq £ M and by vector of initial velocity at the point po. They can be 
used to shift curves in M. Suppose that the function x = x(s) determines parametric 
curve 7 with parameter s in local coordinates x^, . Let's set up a Cauchy problem 
for the equations (1.1) at the points of this curve: 



(1.3) 



lt=0 



x(s), 



lt=0 



z/(s) • n(s). 



Here n(s) is normal vector of the curve, while i^(s) is a scalar function determining 
modulus of initial velocity on the curve 7. The solution of Cauchy problem is a 
vector-function x = x(s,t), it determines the shift of curve 7 along trajectories of 
dynamical system (1.1). Denote it by /t : 7 — > 74. 
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Definition 1.1. Shift ft of curve 7 along trajectories of Newtonian dynamical 
system is called a normal shift if all curves 7t obtained by this shift are perpen- 
dicular to the trajectories of shift. 

This definition doesn't differ from corresponding definition in Chapter II). An- 
other definition from Chapter II), which introduce the concept of dynamical system 
admitting the normal shift, also remains unchanged here. 

Definition 1.2. Newtonian dynamical system (1.1) is called a system admit- 
ting the normal shift in M if for any smooth curve 7 given by vector-function 
X = x(s) and for any point s = sq on 7 one can mark some part of this curve 
containing the point So, and one can find nonzero function i/(s) on this part such 
that shift determined by initial data (1.3) is the normal shift of marked part of 
curve 7 along trajectories of dynamical system (1.1). 

Let's normalize the function i/(s) from (1.3) by the following condition at the 
marked point s = sq: 



Here vq 0. Motivation of such choice of normalizing condition can be found in 
§2 of Chapter II). The condition (1.4) is used in the statement one more concept 
introduced in Chapter II). For present case of dynamical systems in Riemannian 
manifolds it is formulated as follows. 

Definition 1.3. Newtonian dynamical system (1.1) in two-dimensional Rie- 
mannian manifold M satisfies strong normality condition if for any smooth 
parametric curve 7, for any point s = sq on 7, and for any nonzero real number 
j/Q 7^ one can mark some part of this curve containing the point sq, and one can 
find a function i/(s) on this part normalized by the condition (1.4) and such that 
shift determined by initial data (1.3) is the normal shift of marked part of curve 
7 along trajectories of dynamical system (1.1). 

The definitions 1.1, 1.2, 1.3 lie in the base of theory of dynamical systems admit- 
ting the normal shift in Riemannian manifolds. Omitting details of developing this 
theory (see papers [B0I6], [Bol7], and thesis [Shr5]), we formulate its main result for 
two-dimensional case dim M = 2. 

Theorem 1.1. Newtonian dynamical system (1.1) satisfies strong normality con- 
dition in two-dimensional Riemannian manifold M if and only if its force field 
F = F(x, v) satisfies weak normality equations: 



(1.4) 



Ks) 



= 1^0- 



(1.5) 
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2 2 

i=l j=l \ r=l / 

Weak normality equations (1.5) and (1.6) have the same form as the equations 
(1.19) and (1.20) from Chapter III. The operator Vj here acts as differentiation 
with respect to i-th component of velocity vector; this coincides with its definition 
by formula (2.9) from Chapter II: 

dF- 

(1.7) V.ir, = 

The difference of normality equations (1.5) and (1.6) from corresponding equations 
(1.19) and (1.20) in Chapter III consists in definition of covariant derivative Vj: 

(1-8) v.F, = g-Er^^.-EE"^rt«S. 

fe=l fc=l s=l 

Scalar field A, which is defined as scalar product of the vector F with unit vector 
N, is interpreted as projection of vector F to the direction of velocity vector v: 

2 

(1.9) A = Y^^Pi- 

i=l 

Scalar field A from (1.9) defines scalar ansatz 

2 

(1.10) Fk = ANk-\v\J2Pl^iA, 

i=l 

which turns to identity the equation (1.5). Second normality equation (1.6), upon 
applying scalar ansatz (1.10), turns to 

2 / 2 2 

E Vi>l + |v| ^ ^ VgA VrV^A - 
i=l \ q=l r=l 



(1.11) 

^ ' 1 2 \ 

- E ^ V,>1 - |v| E Vi>l Pft = 0. 

r=l r=l / 

This is reduced normality equation. It has the same form as the equation (2.5) 
from Chapter III. However, present case of arbitrary Ricmannian manifold differs 
from Euclidean case M = by formula for covariant derivatives Vj. Therefore 
further simplification of the equation (1.11) require some additional information 
about Riemannian metric in two-dimensional case. 
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§ 2. Isothermal coordinates. 

Next fact is well known in geometry: on any two-dimensional Riemannian man- 
ifold M in some neighborhood of any point p & M there are isothermal local 
coordinates and such that metric tensor g has conformally Euclidean form 



(2.1) 



9io 



-2/ 











-2/ 



in these coordinates. Here / = /{x^jx"^). The proof of this fact can be found in 
[Novl]. Let's write normality equations (1.11) in isothermal coordinates. First of 
all, let's calculate components of metric connection. In order to do this we use 
well-known formula for T^j (see books [Kobl] or [Shr4]): 



(2.2) 



V 9xi 

s=l ^ 



dgis _ dgij 
dx^ dx^ 



Let's substitute components of metric (2.1) into the formula (2.2). This yields 

df , df 



(2.3) 



dxi 



Oik- 



Let's use this formula for to calculate covariant derivatives in the equation (1.11): 

dx^ Z^Z^^^-a^;^' 

fe=i «=i 

2 „ . 2 2 



dx^ dv^ 



Qyk 



yS pfe 



k=l k=l s=l 

Substituting here the expressions (2.3) for connection components, we get 



' dx^ ^ ^ V dx^ ^ dx^ ) dv^ 9v* ^-J ^ dx^ ' 

s=l ^ ' s=l 

dx'^ dv'' dx^ dv^ dx'' dv^ ^-^ dx^ dv^ 
^-^ \ dx^ dx^ I dv^ dv^ dv^ dv^ ^-^ dx^ 

S—l ^ ^ 8—1 



Other derivatives in (1.11) can be calculated without use of connection components. 
We obtain the following formulas for them: 
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Modulus of initial velocity |v| and components of unitary vector N are calciilatecl 
respective to metric (2.1). They depend on conformal factor e"-^, which is present 
in formula for components of this metric: 



|v| = e-f ■ ^(^^TT(^, 



g/ . 



However components of projector operator do not depend on conformal factor 
. They can be calculated in Euclidean metric Qij = Sij: 

Ok (^1)2 + (^2)2- 

Rising one index causes that conformal factor appears again: P'^^ — e^^ P^. Substi- 
tuting all the above expressions for covariant derivatives V,; A, V, Viv4, Vg A, VrVi^, 
and for the quantities A^'', P^, P'^^ into the normality equation (1.11), we get the 
equation for the function A, which contains conformal factor and its derivatives. 
However, if we substitute 



(2.4) A = A'.e-/+^^^rf,.^;'t;^Affe, 

i=l i=l k=l 

then we get the equation with no entry of conformal factor . 

Theorem 2.1. The normality equation (1.11) written for in conformally Eu- 
clidean metric gij = e~^-^ dij is equivalent to analogous equation 

2 / 2 2 

J2 I V^A' + |v| 53 53 P'"^ VgA' VrViA' - 

i=l \ q=l r=l 

2 2 \ 

- 53 iV A' VrV^A' - |v| 53 VrV^A' P^ = 0, 



(2.5) 



r=l 



where covariant derivatives and parameter , P^, P''' are calculated in Euclidean 
metric gij = Sij, while A is expressed through A' according to the formula (2.4). 

Proof of this theorem consists in direct calculations. We shall not give it here. The 
equation (2.5) has exactly the same form as the equation (1.11). This is normality 
equation in Euclidean metric gij = Sij . Its solution A' determines dynamical system 

(2.6) x = v, v = F'(x,v), 

force field of which is calculated by formula of scalar ansatz (1.10) for Euclidean 
metric gij = Sij. The same solution A' of the equation (2.5) determines the solution 
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A of the equation (1.11) bound with A' by formula (1.11). It corresponds to the 
dynamical system written in non-Euclidean metric 

(2.7) x = v, VtV = F(x,v), 

As appears, formulas (2.6) and (2.7) are different forms of writing the differential 
equations of the same dynamical system. This follows from the formula binding 
their force fields: 

2 2 

Coincidence of dynamical systems (2.6) and (2.7) reflects the fact that two con- 
formally equivalent metrics define the same class of dynamical systems admitting 
the normal shift. In two dimensional case this fact is especially important, since 
arbitrary Riemannian metric can be brought to conformally Euclidean form. 

§ 3. Complex form of normality equation. 

Isothermal coordinates determines complex structure on arbitrary two-dimensio- 
nal Riemannian manifold M. Therefore M turns to one-dimensional complex mani- 
fold. Indeed, if and are isothermal local coordinates in M, then z = x^+ix^ is 
complex coordinate. If x^ and x"^ are other isothermal coordinates and z = x^ +ix'^, 
then in the domain of overlapping these two local maps the transition 2 to 5 is deter- 
mined by some holomorphic function z = z{z). This fact is well known in geometry 
(see [Novl]). Let's add another complex w = v^+iv"^ variable to z. This determines 
complex structure in tangent bundle TM. Transfer from complex coordinates z, w 
to other complex coordinates z, w is determined by formulas 

- -/ x - d^i^) 

Z = Z(Z), W = ; W. 

dz 

Riemannian metric (2.1) is hooked to complex structure of M by conformal factor: 
(3.1) g = e-2/(^'^'rf^d^. 

Metric (3.1) determines a class of dynamical systems admitting the normal shift. 
But, according to the theorem 2.1, this class doesn't depend on conformal factor 
e^. Therefore we conclude: each one-dimensional complex manifold M (which in 
complex analysis is called Riemann surface) is canonically bound with some class 
of Newtonian dynamical systems admitting the normal shift. The presence 
(or absence) of metric g is of no importance for this fact. 

Suppose that M is one-dimensional complex manifold. Let's establish direct re- 
lation between class of Newtonian dynamical systems admitting the normal shift in 
M and complex structure of M, without use of metric g. For this purpose we write 
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the equations of Newtonian dynamical system in form (2.6): 

(3.2) x^=v\ x^=v^, 

(3.3) v^=F\ v^ = F^. 

Let's multiply second equation (3.2) by i and add to first one. Then do the same 
with the equations (3.3). As a result we get complex equations of dynamics 

(3.4) z = w, "W = F{z,z,w,w). 

In (3.4) velocity vector is represented by complex number, while force vector F is 
represented by complex fimction F = F{z, z. w, iT>). 

For dynamical system admitting the normal shift the force field is determined 
by scalar function A due to scalar ansatz (1.10). Let's find complex form of the 
relationship (1.10). With this purpose we write it in components: 



(3.5) 



i^.=^7V,-|v|(P,^|4 + P|— ^ 



Components of unitary vector N and components of projector P in (3.5) should be 
calculated in Euclidean metric gij = 6ij : 



7;^ 7;^ 7'^ 



V(t^l)2 + (t;2)2' ^ (^1)2 + (^2)2- 

Let's substitute these expressions into (3.5). Then multiply second relationship (3.5) 
by i and add to first one. As a result we get the expression for F in complex form: 

/o.x r. w f , dA dA\ 

Formula (3.6) is a complex form of scalar ansatz (1.10). It expresses complex function 
F{z, z, w, w) through real function A{z, z, w, w). 

Now let's find complex form of the normality equations (1.11). Note that in (1.11) 
we have two equations that correspond to k = 1 and k = 2. Let's multiply second 
equation by i and add to first one^. As a result we obtain one complex equation. 



^Surely one shouldn't mix index i and i = \/— T. 
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For the sake of brevity in writing this equation we introduce the operators 



Let's complete operators (3.7) by the following two operators: 

Now we can write normality equation (1.11) in complex form; 

D-A ■ {D-D- - Dt)A - \w\ ■ D-A - 

(3.9) 

- A • (D+ - 1)D-A + \w\ ■ DtD-A = 0. 
While formula (3.6) for the force field can be rewritten as 
(3.10) F=^.(1 + D-)A 

Theorem 3.1. Newtonian dynamical system (3.1) satisfies strong normality con- 
dition in one- dimensional complex manifold M if and only if its force field F is deter- 
mined according to formula (3.10) by some real function A = A{z, z, w, w) satisfying 

the normality equation (3.9). 

Theorem 3.1 is simply the other form of theorem 1.1. It refers to the concepts 
of normal shift and strong normality, which are bound with conformal structure of 
onc-dimcnsional complex manifold M, but in contrast to theorem 1.1, this theorem 
doesn't require the existence of some metric on M. 
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